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ABSTRACT 



Context. Thin X-ray filaments are observed in the vicinity of young supernova remnants (SNR) blast waves. Identifying the process 
that creates these filaments would provide direct insight into the particle acceleration occurring within SNR and in particular the 
cosmic ray yield. 

I Aims. We investigate magnetic amplification in the upstream medium of a SNR blast wave through both resonant and non-resonant 

• . regimes of the streaming instability. We attempt to understand more clearly of the diffusive shock acceleration (DSA) efficiency by 

i-C ■ considering various relaxation processes of the magnetic fluctuations in the downstream medium. Multiwavelength radiative signa- 

tures originating in the SNR shock wave are used to test various downstream turbulence relaxation models. 

Methods. Analytical and numerical calculations that couple stochastic differential equation schemes with ID spherical magnetohydro- 
, dynamics simulations are used to investigate, in the context of test particles, turbulence evolution in both the forshock and post-shock 

■ regions. Stochastic second-order Fermi acceleration induced by resonant modes, magnetic field relaxation and amplification, and tur- 
^ ' bulence compression at the shock front are considered to model the multiwavelength filaments produced in SNRs. The y-ray emission 

is also considered in terms of inverse Compton mechanism. 

Results. We confirm the result of Parizot and collaborators that the maximum CR energies should not go well beyond PeV energies 

■ in young SNRs where X-ray filaments are observed. To reproduce observational data, we derive an upper limit to the magnetic field 
^ ' amplitude and so ensure that stochastic particle reacceleration remains inefficient. Considering various magnetic relaxation processes, 

we then infer two necessary conditions to achieve efficient acceleration and X-ray filaments in SNRs: (1) the turbulence must fulfil 
the inequality 2 - p - > 0; where /? is the turbulence spectral index and is the relaxation length energy power-law index; (2) 
the typical relaxation length must be of the order the X-ray rim size. We find that Alvenic/fast magnetosonic mode damping fulfils all 
fSJ ' conditions; while non-linear Kolmogorov damping does not. By confronting previous relaxation processes with observational data, 

we deducte that among our SNR sample, data for the older ones (SN1006 and G347.3-0.5) does not comply with all conditions, which 
means that their X-ray filaments are probably controlled by radiative losses. The younger SNRs, Cassiopeia A, Tycho, and Kepler 
pass all tests and we infer that the downstream magnetic field amplitude is in the range of 200-300 /i Gauss. 



Key words. ISM: supernova remnants - Physical data and processes: Acceleration of particle - Magnetohydrodynamics (MHD) 
Shock waves - Turbulence - Supernova: individuals: Cassiopeia A - Tycho - Kepler - SN1006 - G347.3-0.5 



H 1 . Introduction second in angular size as reported in I Parizot et al' ('2006V Their 

■ - 1 true width, however has to be inferred from deprojection cal- 

Chandra high-angular resolution X-ray observations of young culations by taking into account Ih e curvature o remnant 

supernova remnants (SNR) such as Cassiopeia A, Kepler or Perezhko et al| | 2003aH Ballet||2Q06V This size should depend on 

Tycho, have detected very thin X-ray filaments, which are the magnetic field strength, local fluid properties (the shock ve- 

probably associated wifli the supernova (SN) f orward shock l°"ty ^"^^ compression ratio), and flie relativistic electron diff^u- 

expanding into the interstellar medium (ISM) (Gotthelfet al' sion regime. 

2001; Hwang et al 2002; Rho et al 2002; Uchiyama et al 2003; Measurements of the X-ray rim size inferred a lower limit 

Cassam-Chenai et al. 2004; Ba mba et al I 12005 at iBamba et al to the magnetic field located downstream from the shock front. 

2005b; .Cassam-Chenai et al..,20()7| ). The physical properties of Typical field strengths of two orders of magnitude abo ve the 



these filaments were reviewed by |Vin k & Laming (2003), VinH standar d ISM values Bqc were reported by e.g. , Berezh ko et all 

(2004'),'Weisskopf & Hughes' (2006^Ballet (2006), Pariz ot etall (2003a),'Vink'(2004), Volk et al (2005), Pa rizot et all(i2006l) . and 

(2006), Bamba et al (2006), and Berezhko (2008). These fila- Berezhko (2008). Vink (2004) showed that advective and diffu- 

ments are believed to be produced by synchrotron radiation emit- sive transport also contributes to the filament extension at high 

tedbyTeV electrons. Rim-like filaments are usually of a few arc- energy close to the electron cut-off. The aforementioned con- 
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straints infer a value of the electron spatial diffusion coefficient 
that is a few times higher than the Bohm limit in the downstream 
region from the shockQ. These results support the standard sce- 
nario of diffusive shock acceleration (DSA) in SNRs and im- 
ply that strong magnetic field amplification occurs at the shock 
precursor. However, Chandra observations have been obtained 
for a limited frequency range. Thus, diffusion regimes differing 
from that of the Bohm diffusion c annot be ruled out by these 
sole observations (Mar cowith et alll2006l) . For instance, alterna- 
tive diffusion regimes may affect high energy particle transport 
and modify the way in which the synchrotron spectrum cut- 
off is reconstructed from the ext rapolation of the radio spec- 
trum (IZirakash vili & Aharonianl2007). However , the hard X-ray 
detec tion of SNR RXJ1713-3946.5 by Suzaku dTakahashi et all 
I2OO8 !) supports a cut-off spectrum in agreement with a Bohm- 
like diffusion regime. 

The origin of the magnetic fluctuations sustaining the diffu- 
sive behavior of non-thermal particles remains widely debated. 
One possibility is that the turbulent magnetic field is generated 
by the relativistic particles themselves by means of their stream- 
ing motion ahead of the shock front ( Bell & Lucek 2001) . The 
field amplification has strong implications for the physics of 
cosmic-ray (CR) acceleration at SNR shocks. For instance, a 
calculation including the effect of non-linear turbulence transfer 
concluded that proton acceleration is possible up to the CR spec- 
trum knee at ~ 3 x lO'^eV. This calculation was performed in 
the most extreme shock velocity regimes, particularly for SNRs 
propagating in a hot interstellar medium free of ion-neutral wave 
damping (Ptuskin & Zirakashvili 2003). Bell (2004) discussed a 
non-resonant regime of streaming instability that can generate a 
very strong turbulent magnetic field (and boost the CR maximum 
energ y) readily at the very early st age of the SNR free exp ansion 
phase.lDiamond & Malkovl (12007 ^ and'Pelletie r et all (12006 ) also 
highUghted the importance of determining the saturation level 
of the magnetic fl uctuations, which w as partially discarded in 
the previous work. iPelletier et all (1200 6) demonstrated that both 
resonant and non-resonant regimes of the streaming instability 
have to be considered simultaneously to fix the magnetic field 
spectrum and strength at the shock front. In fast shocks, the non- 
resonant instability dominates the magnetic field generation, the 
level of fluctuation at the shock being found to be similar to 
the value derived by Bell (2004). The resonant instability dom- 
inates in slower shock regimes. The turbulence generated up- 
stream may then relax downstream from the shock front, limi- 
tating of the s patial extent of the non-thermal particle journey 
dPohl et all 120 05). This possibility has not yet been completely 
taken into account in the DSA process and the corresponding 
maximum energy reachable by relativistic particles. This issue is 
investigated in a dedicated section of the present article. We note 
that the problem of the maximum CR energy was addressed by 
IZirakashvih & PtuskinI (^2008) using a semi-analytical approach 
to the non-resonant streaming mode generation. The authors 
identified the maximum CR energy, between the two confine- 
ment limits, expected for a standard ISM medium or a com- 
pletely amplified magnetic field. One should keep in mind that 
several effects can alterate these conclusions s uch as the prop- 
agation into a partially ionised medium (Bvkov & Toptygliinl 
l2005HReville et alll2007l) . thermal effects in the dispe rsion rela- 
tion of the non-resonant instability ("Reville et al'2008*), or a back 
reaction on the CR current ( Riquelme & Spitkovskv 2009) . 



' The Bohm diffusion coefficient is obtained when the particle mean 
free path is equal to its Larmor radius = EjZeB, i.e., DBohm = ''lc/S. 



Al thoug h disputed (see dis cussions in iKatz & WaxmanI 
(2008*) and 'Morli no et all (|2008')), the production of relativis- 
tic hadrons in SNRs is consistent with detection of a few TeV 
y-ray emitting SNRs in the Galactic plane by the HESS tele- 
scope. This y-ray emission may favor the interaction of relativis- 
tic hadrons with a dense molecular c loud leading to the Compton 
upscattering of low energy photons (Aharo nian et alll2bo4l200d 
Al bert e t al 2007). Nevertheless, more observations are manda- 
tory before drawing any firm conclusion about this important 
issue. 

The present article investigates DSA processes involving 
magnetic field amplification and relaxation. The paper considers 
the effect of shock acceleration, spatial variation in the magnetic 
field (and the coiTesponding diffusion coefficient), the possibil- 
ity of finite diffusive extension zones, and the effect of stochastic 
Fermi acceleration by the electromagnetic fluctuations generated 
in the shock precursor. This modelling is performed by means 
of numerical calculations. The numerical scheme is based upon 
the stochastic differential equations (SDE) and is described in 
Appendix|C| Sect.|2|presents the general framework adopted in 
this article. In particular, it investigates the conditions required to 
develop turbulence upstream from the shock, as expected from 
the non-linear evolution of the various regimes of the streaming 
instability. Sect.|3]and|4|investigate the impact of post-shock tur- 
bulence upon particle acceleration. Sect.|3]dealing with advected 
downstream turbulence and sect.|4|refering to a downstream re- 
laxing turbulence. All calculations are the n compared with t hose 
for a sample of young SNRs presented in lParizot et all (l2006l) al- 
ready. 

Table ([TJ summarises the notations used in this article (the sec- 
tion where the parameter is reported at first is also indicated). 



2. Upstream turbulence generation and accelerated 
particle diffusion 

Highly turbulent supernova shocks involve several complex pro- 
cesses that modify the standard DSA model at some stage of 
the SNR evolution. In the upstream region, the properties of the 
turbulence are contr olled by the fastest growing instability and 
its saturation mode (IPelletier et alll2006l) . The diffusion regime 
strongly depends on the competition between the wave growth 
and the e nergy transfer to other scales provoked by non-linear 
cascades (iMarcowith et alll2006l) . The turbulence is then com- 
pressed at the shock-front, i.e., parallel modes (parallel to the 
shock normal) have wavelengths that are shorter by a factor 
equal to the (sub)shock compression ratio. In the downstream 
region, the t urbulence can either be relaxed (Poh l et al 20051) 
or amplified (IPelletier et af 2006: Zir akashviU & Pfii skin 200^. 
The turbulent magnetic field coherence length may also vary 
with the distance t o the shock, whi ch can be modelled using self- 
similar solutions (iKatz et alll2007l) . 

Section 12.11 summarises the properties of the two regimes 
(both resonant and non-resonant) of the streaming instability 
as well as the magnetic field profiles inserted into the coupled 
SDE-magnetohydrodynamics (MHD) numerical calculations. In 
sect. 12.3! we derive the general form of the diffusion coefficient. 
Finally, sect. I2.4l displavs the general expression of the particle 
distribution function, at the shock front, expected in the case of 
spatially varying diffusive zones. The various expressions de- 
rived in this section will be used in sect.|3]and|4| 
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Turbulence parameters 





One D power-law spectral index of the turbulence spectrum [Eq.ll6)] 


m 


Level of magnetic fluctuations with respect to the mean ISM magnetic field [Eq.ll6J] 


<f> Logarithm of the ratio of the maximum momentum to the injection momentum [Eq.(|4J] 


^max 


Longest wavelength of the magnetic turbulence spectrum (Sect. 12. 3b 




Coherence length of the magnetic fluctuations (Sect. 12. 3t 


CT 


Normalisation factor entering the turbulent spectrum (Sect. 12. 3b 


Saji 


Power-law energy dependance index of the relaxation lengths either up- 
or downstream (Sect.|4]and [Eq.ll23l) 


H 


Ratio of the upstream to the downstream diffusion coefficient at the shock front [Eq.(ll6b] 


Sb 


Ratio of the resonant to the non-resonant magnetic field strength at the shock front [Eq.l|3J] 



Relativistic particle 
parameters 





Ratio of the CR pressure to the shock dynamical pressure [Eq.(|3J(] 


rt 


Larmor radius of a particle (defined using resonant magnetic field) 


P 


Ratio of the particle Larmor radius to /i„,ax/27r (also called reduced rigidity, see sect. 12. 3t 


EcR-mux 


Maximal cosmic ray energy (Sect. 12.311 


p 


Maximal electron energy (Sect. 13. U 




Cut-ofi' synchrotron photon energy emitted by electrons at iie-max (Sect.|3JJ 


£cR-min 


Injection energv of the cosmic ravs (Sect. 14. 1.4b 




Energv of the electrons producing the observed X-rav filaments (Sect. 14.2.4b 



SNRs parameters 



Kh,4 


Velocity of the SNR shock wave (in 10'' km/s unit) 


5d/u,-4 


Magnetic field amplitude at the shock front respectively in the 
down- and upstream medium (in lO"'* Gauss unit) 


''Bj ''sub I ''tot 


Magnetic, sub-shock, and total shock compression ratios (Sect.l3.1b 


Ai?x.-2 


X-rav filament deprojected width (in 10"*^ parsec unit. sect.l4.2.4b 



Equation parameters 





3r'/(r- l)[Eq.(19)] 




q(fi) X (H(r,p)/r + I) [Eq.(36)] 


f 

^sync 


H(r,p) + r/H(r,/3)/r'^ + r [Eq.(39)] 


g(r) 


3/(r- l)x(//(r,yS)/r+ l)[Eq.(40)] 


C(6d) 


(£e-ntax/£.-«fo)''" [Eq.(41)] 



Table 1. Summary of the notations used in this article to denote the various physical quantities and parameters involved in our 
description of high energy particle yield in supernova remnants (SNR). 



2.1. Cosmic-ray streaming instabilities 



The streaming instability which is provoked by the superalfvenic 
motion of accelerated energetic particles, generates magnetic 
fluctuations over a large interval of wave numbers. The res- 
onant instability involves wave-particle intera ction on wave 
scales of the order of the particle gyro-radius (ISkillingll 19751: 
iBelT & Luce k 200 il).. The non-resonant re gime was adapted t o 
SNR shock waves by Bell (2004) (see also Pelletier et al ( 2006|), 
IZirakashvili"& Ptuskin (2008) . and Amato & Blast (2009.) for 
further details). The non-resonant waves are produced, at least 
in the linear growth phase of the instability, on scales much 
smaller than tl. However, the ability of the instability to both 
enter deeply into the non-linear regime and saturate at a mag- 
netic field level 6B » Boa remains debated ( Reville et"all2008t 
iNiemiec et all2008t iRiquelme & Spitkovskvll2009l) . In the next 
paragraph, we summarise the main properties of the wave 
modes generated by the non-resonant streaming instability (sect. 
( 12.1.1b ). We then present the characteristics of the resonant 
regime in sect. ( I2.1.2l i. 



2.1.1. The non-resonant regime 

In the linear p hase, the m ost rapidly growing waves have large 
wave numbers (lBellll2004l) given by 



k <kr 



(1) 



where j„ - ncr^Vsh is the current produced by the cosmic rays 
ahead of the shock wave, ncR is the CR density, and Vsh is the 
shock velocity measured in the upstream restframe. 

The wave number corresponding to the maximal growth rate 

Tmax 



■ ^up^aoo is 



^up — ~T — X 

2 rico 



■cp ■ 



2V^ 



(2) 



where Woo, f^cp = eBoo/('«pc), and Vaco = Boo/ yAmiZ are, re- 
spectively, the density, the cyclotron frequency, and the Alfven 

velocity in the ISMB 

MHD calculations (lBellll2004t IZirakashvili & Ptuskinll2008h 
have shown that beyond an exponential growth phase located on 
typical scale of 



Vsh/r( 



CR-max 



^ The density n^o is usually the ion density, but when the coupling 
between ion and neutrals is effective it must also involve the density of 
neutrals 
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from the shock, the instability enters the non-Hnear regime. 
The magnetic fluctuations are redistributed on larger scales, 
while the turbulence level evo lves in a linear way. iBelll 
(12004 and iPeUetier eta 1 d2006h discussed several saturation 
processes that all lead to an energy transfer from the dom- 
in ant wavelength towards long wavelengths (see discussion 
in iRiquelme & Spitkovskvl feOOQ)). One may then expect the 
coherence length of the turbulence to be transferred from a 
scale 4oh-L, where k^^^ < 4oh-L <K '^L-CR-max, to a scale 

4oh-NL < ''L-CR-max whcrC TL-CR-max - '"L-CR-max X Boo/B is 

the renormalised maximum energy CR gyro-radius in the ampli- 
fied magnetic field B. Resonant modes have a harder spectrum 
(iPelletier et al 2006) hence contribute to an increase in the co- 
herence length of the turbulence (see sect. 12.1.21 ). So, hereafter, 
we consider both regimes producing a turbulence with a coher- 
ence scale close to pL-cR-max, i-e., we neglect the extension of 
the upstream region where the non-resonant instability is in the 
linear regime (see sect. 12.31) . 

Another important property of non-resonant mode s is that 
they have non-vanishing helicity (IPelletier et all 120061) . These 
modes are mostly proton-induced and have a right-handed po- 
larisation with respect to the mean magnetic field far upstream. 
This non-zero helicity may be the origin of additional amplifica- 
tion in the downstream medium, where the total magnetic field 
can eventually reach values close to the equipartition with the 
kinetic energy of the thermal gas. 



2.1.2. The resonant regime 

The resonant re gime develops simu ltaneously with the non- 
resonant regime ( Pelletier et alll200"6l) and cannot be discarded. 
The total amplification factor of the magnetic field A^^^ - 
^tot/^i a distance x from the shock front is a combinaition of 
both non-resonant and resonant contributions, namely A^^^ix) - 
A^^(x)+A^{x). The exact spatial dependence of Ar(x) is derived 
in Appendix|A|for completeness. It is found that a good approx- 

1 /2 

imation is Ar oc A^^. 

To quantify the previous assert ion, we parametrise t he contribu- 
tion of each instability regime. Pelletier et a l (l2006l) argued that 
the shock velocity is the main controlling factor of each con- 
tribution. This dependence can be inferred from Eq.( IA.2b . By 
comparing the respective saturation values of each regime, one 
finds that 

Br(x^O) ^ ^ I^CRcV'" 
Bnr(x = 0) I / ' 



(3) 



while 



Bnr{x = 0) 

Bca 



CR 



1/2 



(4) 



The level of magnetic fluctuations at the shock front given by Eq. 
(O and (HJi is controlled by both Sb and the fraction ^cr of the 
SNR dynamical pressure transferred into the CR. The parameter 
(f) - log(/:imax/ Pinj) is the logarithm of the ratio of the maximum 
to the injection momentum and is approximately between 15 and 
16. 

As a fiducial example, we assume that ^cr - 0.2, B^o - 
4 yuGauss, and that the ion density asni - 0.7 cm"^. We then 
identify three distinct domains: 

1. (5b > 3 (corresponding to V^h ^ c/400) in which the mag- 
netic field amplification provided by the streaming instability 
is modest for slow shock velocities. 



2. 1 < 5b < 3 (corresponding to c/400 < Vsh < c/10): for 
which we get the ordering Br > Bnr > ^oo and, that the 
ratio Br/Bnr does not exceed a factor 2. 

3. 6b < I (corresponding to Vsh > c/10): for which the mag- 
netic ordering becomes Z?nr > Br > Bca- In that case, an 
upper limit velocity stands close to c. Beyond that limit, the 
amplification by the non-resonant instability is maximal. An 
accurate analysis is then necessary to compare the saturation 
of the i nstability induced b y both advection and non-linear 
efl^ects (IPelletier et alll2009l) . 

Electrons and protons (or ions) moving in the forward or 
backward direction can resonate with either forward or back- 
ward modes. Efficient mode redistribution is expected to pro- 
duce waves i n both directions in the shock precursor (see the 
Appendix of IPelletier et all (l2006l) for further details). We note 
that the interaction between resonant Alfven waves and the 
shock produce magnetic helicity that is different from either 
+ 1 or -1 and ensures that second order Fermi acceleration by 
the resonant modes is unavoidable in the downstream regio n 
dCampeanu & Schlickeiser 1 [T99irV^inio & Schlickeise3[T999h . 
This effect is discussed in sect. 13.21 



2.2. A note on the evolution of non-resonant modes 

Non-resonant modes are purely growing modes of null fre- 
quency, at least in the linear phase. They do not correspond to 
any normal mode of the plasma as in the case of the resonant 
regime. Consequently they are expected to be rapidly damped 
once the source term is quenched, i.e., at the shock front. The 
damping length should be of the order of a few plasma skin 
depths. However, these modes also have a non-va nishing helicity 
(IPelletier et alll2006l: IZirakashvili & Ptuskinll200 8) (as we see in 
sect. 14.1.41 ). So a fraction of the turbulent spectrum can grow fur- 
ther downstream by means of dynamo action. At this point, the 
downstream evolution of the non-resonant spectrum is unclear. 
In some conditions the combinaition of magnetic field compres- 
sion and non-resonant mode damping at the shock front leads to 
a downstream magnetic field that is weaker than the upstream 
field, especially in the very fast shock regime (regime 3. dis- 
cussed in sect. 12.1.2b . This is not the case for the SNR sample 
considered in this work as the resonant modes tend to be dom- 
inant at the shock front. A complete investigation of this diffi- 
cult issue would require a detailed investigation of the interstel- 
lar medium interaction with shocks to fix the ratio Br/B^r. For 
this reason, we assume hereafter that the downstream behaviour 
of the turbulence is dominated by the resonant mode. However, 
even if Br/Z?nr > 1 at the shock front, the fastest growing chan- 
nel is the non-resonant one, which is important for the complete 
setting of the magnetic field turbulence in the upstream region. 
We acknowledge that this assumption weakens the analysis pre- 
sented in the following sections and consider this first work to be 
an attempt to isolate the main properties of the turbulence around 
a SNR shock. 



2.3. Upstream diffusion regimes 

As previously discussed, the most energetic CRs generate fluc- 
tuations at scales that are much smaller than rL-cR-max- These 
particles experience small-scale turbulence exclusively in the un- 
amplified magnetic field. Thus, the diffusion coefficient at max- 
imum energy scales as D(£'cR-max) - (^l -CR-max / 4oh) ^coh^ 
(see below). This allows us to compare Xg and ^ditf(£^CR-max) = 
D(£'cR-max)/Vsh, the diffusivc length of the most energetic cos- 
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mic rays. One can then write (iPelletier et alll2006l) 



Vsh 
X X 



^coh 



2(f> Pboo 
3^CR PcoVl. Vaoo ^L-CR-n 



X ^diff(£'cR-max) ■ (5) 



We find that Xg <k ^diffC^CR-max) in fast shocks (Vsh > 10 ^c) 
because (Pboo/PooV^^) Vsh/Vac« ~ Vaoo/Vsh- The following no- 
tations were used to derive the previous result: the CR den- 
sity is linked to the CR pressure by ncR - ^Pcr/<Pp*c and 
P* - fCR-max at a distance x - ^diff(£^CR-max) from the shock. 
The parameter ^CR is probably between 0.1 and 0.3. 

CRs and electrons of energy lower than EcR-max, diffuse by 
means of a large-scale turbulence, their t ransport properties dif- 
ferin g from those of most energetic CRs jZirakashvili & PtuskinI 
12008 ). Whatever the turbulence level, the angular diffusion fre- 
quency (for a re lativistic pa r ticle i n an amplified field) can be 
estimated as (see lCasse et all ( l2002h . their Eq.A22) : 



-r^^x{j3-l)xbc 



6B^ 
W 



(6) 



where 



fe = 4ohX flflnW(Mcoh)"^. (7) 

•^^'miii-NR^oh 



The turbulence spectrum is assumed to spread over the range 
A;~ jjj] with a ID power-law spectral index fi.lf fi - 1, the 
term 1/06 - 1) has to be replaced by a factor cr - ln(A;max/^min)- 
The corresponding spatial diffusion coefficient is by definition 
Z) = c^/3vs. Its energy dependence is related to the development 
of the instability. In the linear phase (small scale turbulence), we 
recover the above expression for ^diif(£cR-max)- If ^min-NR^oh - 
1, after introducing the level of turbulence 777 - SB^/B^, we find 
that 

DiE)^-f-.'-^.ljl-]\ (8) 
n(J3-l) T/T \4oh/ 

In the non-linear phase (i.e., large-scale turbulence), we have 
^min-NR ~ ''l and so 



D(E) 



P 4ohC 

^ ^ ^ 

n(Ji-\) jjT \4oh 



(9) 



The results obtained bv lCasse et all ( |2002|) can be recovered us- 
ing {coh = pM^m-dx/^n by adopting a reduced rigidity of p = 
27rrL(£')//lmax- The length /i^ax - 204oh is the maximum scale 
of the turbulence and pM is a number ~ 0.2 - 0.3. This latter 
number corresponds to the reduced rigidity at which the transi- 
tion between the two diffusive regimes operates. For instance, 
assuming 77T - 1 and /3 - 5/3, one finds that D ^ 2.2DBohm 
at Pl = ^coh, which i s cons istent with the numerical solutions 
found by ICasse et a 1 (I2OO2I) . IfyS = 1, the energy independent 
ratio D/DBohm - 3o-/n ^ 15 - 16. 

We hereafter refer to q(p) as the normalization of the diffu- 
sion coefficient such that 



Z)(£).^X^^x(ii- 



l-f} 



It is noteworthy that the normalization of the diffusion coeffi- 
cient is given by qifi) and must not be confused with the nor- 
malization of the turbulent spectrum. Both quantities appear to 
have similar expressions as seen from quasi-linear theory cal- 
culatio ns or from numerical estimates obtained in ICasse et all 
(I2OO2I) . Nevertheless, they differ in a strong turbulence regime. 



iRevilleet all(l2008h discussed some solutions by clearly display- 
ing diffusion coefficient with sub-Bohm values. This issue is be- 
yond the scope of the calculations performed here and its con- 
sideration is_postponed to a future work (see also a recent work 
bv lShalchl (12009.) ). Considering these uncertainties, we consider 
q(l3 ) as a free parameter hereafter. 

IPelletier et all (|2006) obtained a ID stationary P - 2 power- 
law solution for the non-resonant wave spectrum. We can see 
from the above analysis that the energetic particle transport prop- 
erties around the shock front depend on the possibility that non- 
resonant instability will deeply enter in the non-linear regime. 
Verifying this condition leads to a diffusion coefficient at £ <s; 
£^CR-max givcn by Eq.(|9]l,the magnetic field profile being char- 
acterised by an exponential growth over a scale Xg and a linear 
growth over a scale x < ^diff(£^CR-max)- 

This qualitative analysis confirms that the non-resonant in- 
stability contributes to the turbulence level over a wide range of 
parameters (once the non-linear regime of the instability is es- 
tablished) and the control of the tur bulence coh e rence length. 
However, the analysis presented in P elletier et all (|2006') shows 
however that the resonant instability at least in the domain 2 
of our fiducial example above also contributes to the magnetic 
fluctuation spectrum. The resonant wave spectrum is found to be 
harder, i.e., for a CR distribution spectrum scaling as p"^, the ID 
turbulence spectrum has an index fi - 1 . In this work we assume 
that the turbulence index is in the range 1 < /3 <2. 

2.4. Shock particle distribution 

Before discussing the effect of turbulence evolution in the down- 
stream region, we present the general solution for the particle 
distribution at the shock front in the case of spatially varying 
diffusion coefficients, where radiative losses are discarded. The 
complete calculation is presented in AppendixiBl We briefly out- 
line our result (see Eq. IB. 51 ) as follows. We have assumed up- 
stream and downstream magnetic fluctuations variation lengths 
^u/d to be scale (or energy) dependent (see sectionlDi. The slope 
of the stationary particle distribution (neglecting any radiative 
loss) at the shock front is 



dlnfsip) 
dlnp 



3r 



D,(0,p)exp(^C^e,(x',p)dx'^ 
Uu J^^ exp(J^'^ 6a(x' , p)dx'^ dx 
Z)d(0, p) exp ^- J^' 0u(x', p)dx'^ 



rud J^'' exp ^- J^'' 9d(x' , p)dx'^ dx 



(11) 



The value of the spectrum slope is controlled by the functions 
6'u/d - Mu/d/i'u/d - dlnDu/d/dx (see Eq lB.3l l. In the basic case 
where both upstream and downstream diffusion coefficients can 
be assumed as space independent over lengths /"u/d from the 
shock (and to v anish beyond these distances), th e above expres- 
sion reduces to (lOstrowski & Schlickeiseilll996l) : 



(10) dlnfs(p) 



din p 



1 



r - 1 \ 1 - exp(-Mu4/£>u) exp(Md4/-Dd) - 1 / 

(12) 

If the shock wave is modified by the CR back-reaction, r will 
depend on the particle energy and the shock spectrum will not 
behave like power law. We note that provided functions 6 are 
large compared to unity, the previous relation indicates that we 
obtain standard power-law spectrum expected from DSA theory. 
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The present article investigates the effects of energy and spa- 
tial dependencies of the 6 functions in both up- and downstream 
regions, by relying on a set of available multiwavelength data of 
five SNR: Cassiopeia A, Tycho, Kepler, SN1006, and G347.3- 
0.5 (also known as RXJ 1713-3946.5). All of these remnants are 
of the case 2 discussed in sect. 12.1.21 and correspond to mildly 
fast shocks where both resonant and non-resonant magnetic field 
amplification occur. 

3. Particle acceleration in the case of a downstream 
advected magnetic field 

This section examines the DSA process for an efficient tur- 
bulence amplification mechanism producing a strong magnetic 
field in the shock precursor (see sect. |2|. In the first sect. ( 13. lb . 
we reconsider the calculations performed by Parizot et al (200q) 
but this time including the effect of turbulent scale compression 
at the shock front. Section ( I3.2l i then addresses the usually over- 
looked aspect of stochastic particle acceleration in the down- 
stream flow. Finally, sect. ( I3.3l l deals with tests involving the 
shock solutions obtained by Zirakashvili & Aharonian ( 2003) 
for various turbulent spectrum scalings. We then incorporate par- 
ticle losses and Fermi stochastic acceleration into the Fermi cy- 
cles and proceed with different numerical experiments. We con- 
clude with a comparison between X-ray and y-ray filaments pro- 
duced by inverse Compton up-scattering of cosmic microwave 
background photons. 

3. 1 . Downstream diffusion regimes and maximum particle 
energies 

Downstream of the shock, the particle distribution was fully 
isotropised (to an order of V/v) and the streaming instability 
quenched. We insert the magnetic profiles derived in the pre- 
vious section into the diffusion coefficients (see Eq|9]l. To de- 
rive the downstream diffusion coefficients, we need to spec- 
ify properly how the transition occurs at the shock front. We 
only consider here the case of a strong magnetic field amplifica- 
tion at the shock precursor. The upstream magnetic field being 
highly disordered, the magnetic compression ratio then becomes 

^(l+2r2j/3 < r,ub (with r.^b > ifl 

-1/2 



Bu = fid X 



1 +2 



sub 



Bd 



(13) 



iParizot etal (l2006l) only considered this final effect. But in the 
meantime, the maximum turbulence scale downstream is re- 
duced by a factor r^^h: 

'tmax-d - ■ {^^) 

/"sub 



^ We make a distinction between the compression ratio at the sub- 
shock (^siib < 4) and the total shock compression ratio rtot > 4. In the 
case of weakly modified shocks, we have rt„t - ^sub - = 4. In the 
case of strongly CR modified shocks, one obtains rtoc > r > r^ub- If the 
sole adiabatic heating of the precurso r is considered, value s '^iib = 2-3 
and rtot > 10 are possible (see e.g.. iBerezhko & EllisonI ( 1 19991) ). If a 
substantial gas heating in the precursor is produced for instance by the 
absorption of Alfven waves, the total compression ratio ca nnot be muc h 
higher than 10, under ISM conditions considered above (Bykov 2004"). 
In a strongly modified shock, the most energetic electrons producing 
the X-ray filaments have energy E » EcRmin and do experience a com- 
pression ratio close to rtot. This value will be used in the following es- 
timations. Values of r^^f, = 2 and rf„t = 10 are accepted in this work in 
the case of strongly CR modified shock. 



This scale compression induces an enhancement of the tangen- 
tial magnetic field component and a reduction in the maximum 
turbulence length in the downstream region. The downstream 
turbulence is then anisotropic, display ing elongated e ddies in the 
direction parallel to the shock front (Marcowith et al 2006) un- 
less other non-linear processes prevail (Zirakashvili & Ptuski^ 
2008). The coherence length of the turbulence is hereafter as- 
sumed to be a constant. 

We can define the downstream diffusion coefficient accord- 
ing to the definition of the upstream coefficient given in Eq.(|9l) 



— X X I 

n 27r77T-d \Pm 



(15) 



In the remaider of the article, we only consider the case where 
j8u=y3d=y6. 

Using Eq.(|9l) evaluated atx - Q as well as Eqs. ( flj] ) and (fT4l i. 
we end up linking up- and downstream diffusion coefficients at 
the shock front (where we have assumed that rij ^ 1): 



i ^ sub 



2-/3 



= Dd X //(rsub,y8) , 



(16) 



Once the up- and downstream diffusion coefficients are set, the 
magnetic field at the shock front c an be inferred follo wing the 
same procedure as that adopted in iParizot etal ( l2006h (see the 
article for the detailed derivation). The balance between the elec- 
tron acceleration rate and the mean synchrotron loss rate fixes 
the maximum electron energy to facc(£'e-max) - (fsyn(£'e-max))- 
The synchrotron l oss timescale is obtained from Eq.(17) of 
IParizot etall (l2006b using the mean square magnetic field expe- 
rienced by relativistic electrons during one Fermi cycle: 



<B'> = Bd X 



H(J3) + r,o, 



Following DSA standard theory the acceleration rate is 



^acc(^) — 



3r2 Dd{E) 



r- 1 



H(r,/3) 



HE) 



+ 1 



(17) 



(18) 



Basic analytical relations can be derived when Bohm diffusion 
regime conditions prevail. In that case, electron and proton ac- 
celerations are no longer related because the diffusion coeffi- 
cient no longer de pends on /l,nax anymor^ Equation (30) in 
IParizot etall (l2006l) can be used to derive the downstream mag- 
netic field amplitude and an estimate of the synchrotron photon 
energy cut-off 



■^y-cut 



[0.875 keV] x 



'^sh,4 



(19) 



where we note that yir^od = ^rljirioi - 1), Aot = not/4, and 
q - qip - 1)/16. The maximum electron energy is found to be 
approximately 10 TeV in our SNR sample, a value close to the 
maximum CR energy. To derive such result, we assumed that the 
compression ratio r at £ = £e-max is approximatively ~ rtot. 

In Table (|2]l, we list the inferred values of the downstream 
magnetic field in the context of an advection dominated X- 
ray rim, where a Bohm-type turbulence is occurring. We have 
also provide the theoretical values of fiy-cut required to verify 
facc(£'e-iuail p_i£awi(£'e-max)>- The similar to those in table 1 of 
Parizot et all (l2006l) . except for SN1006 where we used the value 



excepted at the highest energies. 
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Supernova remnant 


Bi (mG) 


^y— cut 
'-'y— cut,obs 


5d-Fii (inG) 


Cas A 


558 


0.2 


2.7 


Kepler 


433 


0.3 


2.3 


Tycho 


586 


0.7 


1.5 


SN 1006 


170 


0.07 


0.56 


G347.3-0.5 


131 


0.05 


2.1 



Table 2. Inferred values of the downstream magnetic field amplitude and synchrotron photon cut-off energy in the case of an 
advection-dominated rim where Bohm diffusion regime prevails (JS - I and qifi = 1) = cr). The magnetic field values were 
calculated by assuming rtot = fsab - 4. 



of shock velocity (4900A:m/s) given in lAcero et all ( l2007b . The 
results presented in this table were performed using a diffusion 
coefficient normalization qifi - \) - cr corresponding to predic- 
tions by the quasi-linear theory. 

Based on the aforementioned assumptions, it appears that 
older SNRs (Tsnr > 1000 yr) should have a synchrotron cut- 
off energy that is much lower than the observed value. However, 
as for instance in the case of SN1006, the cut-off frequency de- 
pends on the observed region of the SNR and 3 keV is probably 
an upper limit. On the other hand, young SNRs (Tsnr < SOOyr) 
exhibit, in the same context, strong magnetic fields and syn- 
chrotron energies cut-off close to the cut-off deduced from the 
observa tions. The effect is even stronger in the case of modified 
shocks. IParizot et al (2006) noted that the Bohm regime does 
not allow the DSA theory to reproduce accurately the X-ray fil- 
aments unless the diffusion coefficient normalization is replaced 
by a factor of a few. This is confirmed by the close agreement 
between the two cut-off energies obtained for the young SNR. 

Several uncertainties may shift the value of the cut- 
off frequency from the extrapo lation using the radio data. 
IZirakashvili & AharonianI (l2007h pointed out that the electron 
particle distribution can be cut off in a smoother way than by a 
pure exponential cut-off. In that case, the true cut-off frequency 
is shifted towards higher energies. In the meantime, the observed 
synchrotron cut-off used previously is probably to be an upper 
limit because of the back-reaction of CR on the shock structure 
producing a curved shape of the spectrum. It seems justified to 
develop a detailed non-linear calculation to improve the estimate 
of the discrepancy between these solutions with a simple expo- 
nential cut-off. This aspect should also be an important issue for 
the next hard X-ray satellites generation such as nuStar or Next. 
We postpone its investigation to future work. 

To summarise, we can say that scale compression has a very 
limited impact on th e abov e calculation and that the results de- 
rived in lParizot et all (l2006h are quite robust. 

3.2. Considering downstream stocliastic Fermi acceleration 

The downstream magnetic field amplitudes derived in sect.|3]are 
lower limits, while the observed filament sizes are just upper 
limits because of the lack of resolution of X-ray instruments. If 
the downstream magnetic field reaches values close to mGauss 
and does not relax rapidly, then at some stage the Alfven ve- 
locity will be of the order of the downstream fluid velocity. In 
that case, stochastic Fermi acceleration can no longer be ne- 
glected. Electrons will interact with turbulence modes generated 
by the resonant streaming instability since non-resonant modes 
are right-handed polarized and thus cannot interact with elec- 
trons. We included in our numerical calculations the so-called 
Fermi second-order process (in addition to the usual first-order 



acceleration) combined with energy losses, namely synchrotron 
losses for the electrons. We implicitly assume in the remainder 
of the paper that an efficient redistribution of forward and back- 
ward waves operates by means of non-U near interaction with 
magneto-sonic waves (Pelletier et af 2006"). In that case, forward 
and backward modes transmitted downstream are in balance 
dVainio & Schhckeiseilll999i) . This assumption enables us to es- 
timate the magnetic field amplitude when dominant stochastic 
Fermi acceleration occurs. Issues dealing with imbalanced mag- 
netic turbulence are beyond the scope of this paper and will be 
investigated in future work. 

The acceleration timescale characterising the stochastic Fermi 
process for a relativistic particle can be written as 

9D(E) 

faccFII - —775 ■ (20) 

The conditions in which a stochastic acceleration less efficient 
than the usual shock acceleration can be transposed into a condi- 
tion on the downstream magnetic field by writing facc.Fii ^ facc.Fi- 
Using Eqs. ( fTSl ) and ( l20l i. one can readily find that 

1/2 -1/2 

Bd-Fii < [714 /.Gauss] x _, ^. . , (21) 

In this expression, we exceptionally use a shock velocity ex- 
pressed in units of 10^ km/s and ISM density in units of 
0.1 cm-\ 

In the case of young SNRs propagating into a standard ISM 
medium with typical hydrogen densities ~ 10 ' cm"^, the previ- 
ous limit leads to magnetic field strengths ~ 1 - 2 mGauss for a 
typical shock velocity of the order of 5 x 10^ km/s. This is con- 
firmed by the values of the limited magnetic field strengths given 
tab. |2|for each SNR. The surrounding gas density in most cases 
provides only a crude estimate or is derived from averaged val- 
ue s over the entire remn ants. We have used f or Cas A: = 
1 (i Berezhko eTal l2003 b). Kepler: «oc ^ 0.7 dAharonian et all 
20081), Tycho : rip, ^ .4 (Hughesl lIOOOh . SN1006: Hoc = 0.05 
(SB rims see Acero et al (2007)), G347. 3-0.5: «oo = 1 (poorly 
constrained see Aharonian et al (2006 )). 

The Fermi stochastic acceleration process produces an en- 
ergy gain in the downstream medium and a hardening of the par- 
ticle distribution at the shock front (see Eq. 15 in Marcowith et al 
(l2006i) and the simulations in sect. l3.3.2l i. As particles are contin- 
uously reaccelerated downstream, they are expected to produce 
larger X-ray filaments. Both effects seem clearly incompatible 
with the available data. The magnetic field fluctuations in reso- 
nance with electrons are then expected to saturate at the shock 
front with magnetic field amplitude <s; Bd-Fii, which is below the 
value for equipartition with thermal pressure of the flow. 
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3.3. Numerical experiments 



The SDE method presented in Appendix |C] does not account 
for the back-reaction of CR over the fluid flow. This would re- 
quire a special smoothing and a difficult treatment of the CR 
pressure Pgr- The latter calculated from the particle distribu- 
tion f{p, r) at each grid point would produce unphysical fluc- 
tuations that develop with time. Several numerical works have 
started to included wave gene ration effects i n CR modified 
shock hydrodynamics ( Vladim irov et all l2006t iKang & JonesI 
l2007t IVladimirov e tal 2008). Some semi-analytical works has 
also started to investigate the effect of the wave precurso r heat- 
ing on the CR back-reaction process (ICaprioli et aJ2008al) . Both 
approaches seem to reach a similar conclusion: the heating of 
the precursor by the wave damping reduces the gas compress - 
ibility and thus reduces the shock compression (Bvkov 2004|). 
Stationary solutions are found to be rather close to the test par- 
ticle case. Calculations performed in the test particle framework 
using SDEs can then reproduce the main properties of the par- 
ticle acceleration process. SDE have several advantages: they 
are simple to implement and rather simple to couple with MHD 
equations. SDE schemes enable a fast and large investigation of 
the parameter space of the DSA mechanism. For instance, the in- 
clusion of Fermi stochastic acceleration is rather simple in both 
the SDE scheme and in the use of various spatial diffusion coef- 
ficient regimes. Our results can, for instance, be used as limiting 
tests for future non-linear simulations. 



3.3.1. Synchrotron spectrum solutions 



We first validate the aforementioned numerical scheme by 
achieving calculations in different config urations, such as re- 
prod ucing the analytical results of Zir akashvili & AharonianI 
( 1200 7). In this work, the authors define the relativistic electron 
energy spectra at the shock front in the presence of a discontinu- 
ous magnetic field (the discontinuity is located at the shock). We 
performed several SDE-MHD simulations where constant up- 
stream and downstream magnetic fields prevail (Z?j/Z?„ - rs - 
VTT, Tsub is set to 4) and where the shock velocity of the flow 
is set to 3000 km/s. The various presented simulations differ 
only in terms of their implemented spatial diffusion coefficients, 
where D - £>Bohm(£^inf) (£/£inf)'^° (the particles ar e injec ted at 
energy isinj = 5 TeV). IZirakashvili & AharonianI (l2007l) pro- 
vided the shape of the electron energy spectra at the shock front 
beyond the energy cut-off fie-max induced by synchrotron losses, 
namely A^(£') oc exp(-(£'/£'e-max)'^"°)- Figure [1] displays three 
simulations with cd - 1 (Bohm diffusion), q-d = 1 /2 (Kraichnan 
turbulence), and 0-0 = (constant coefficient). The result of the 
numerical calculations are displayed using items, while analyt- 
ical solutions of Zirak ashvili & Ah aronian (2007) are displayed 
using solid lines. In the figure the following parameters have 
been used: velocity of the downstream fluid 3000 km/s, com- 
pression factor Tsub = 4) and uniform upstream and downstream 
magnetic field are set {rg — VTT. We have set various diffu- 
sion regime {D oc E''°) while using our new numerical SDE 
scheme described in appendix |C] The agreement between nu- 
merical calculations and analytical profile is good and proves 
that the skew SDE numerical scheme is valid for all kinds of dif- 
fusion regimes and can handle magnetic discontinuities properly 
(see sect. ( IC.2.11 ) for further details). 



Co 
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D= DJe, 



e/e. 



V,,„=3000 
_B, = 300fiG 

ByB,=(ir" 
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Fig. 1. Shock front energy spectra of relativistic electrons pro- 
vided by multi-scale simulations where the MHD part of the 
simulation mimics the behavior of a SNR blast wave. 

3.3.2. Shock particle distribution and second order Fermi 
process 

Figure |2] and [3] show the shock particle distribution and syn- 
chrotron spectra for the parameters corresponding to the con- 
ditions that prevail in the Kepler and G347.3-0.5 SNRs, respec- 
tively. In the case of the Kepler SNR, we use the parameters 
Vsh = 5.4 X 10^ km/s, = 433 ^G and, y6 = 1. Upstream 
density is 0.7cm"^ (Berezhko et al (2006) estimated the density 
to be iioo < 0.7cm"-'). In the case of G347-0.5, we set param- 
eters to be Vsh = 4000 km/s, Ba = 131 uG, /6 = 1- T he aver- 
aged upstream density is 1 cm"^ (lAharonian et alll2b06l) . In both 
cases, the magnetic profiles used in the simulations are also pre- 
sented. The dashed-lin e shows the stationary solution found in 
iMarcowith et all (l2006l) . which includes particle re-acceleration 
in the Fermi cycle. In the upper right panel, the acceleration 
(with the sole regular Fermi acceleration), and both the diffusive 
and downstream residence timescales are displayed. Diamonds 
are obtained using a numerical calculation of the acceleration 
timescale. The slight excess is produced by the stochastic Fermi 
acceleration process. We also display the synchrotron spectrum 
and the magnetic profile around the shock front at f = 4-QQyr. 

The maximum CR energy (and the aspect ratio A:n,ax/^min) 
corresponds to the maximum CR energy limited by either par- 
ticle escapes in the upstream medium or the SNR age limit. At 
EcR-m-dx, the maximal upstream diffusion coefficient allowed by 
the escaping limit is: 



£>(£'cR-max) = ;r X ^sh V^h 



(22) 



The factor x is usually not accuratly defined. An accurate de- 
termination of this parameter would require to be performed by 
non-linear simulations of DSA that include the effect of the tur- 
bulence generation back-reaction on the flow. A fraction of few 
tenth of percent of the SNR radius is usually assumed in theoret- 
i cal calcul ations and seems to be reasonable (Berezhko 199^ 
iCaprioli et alll2008bl) . The normalization ;f = x/^-^ is then ac- 
ceptable in this text. 

It can be seen from Figs|2] and [3] that stochastic accelera- 
tion slightly modifies the shock particle spectrum in the case of 
the Kepler SNR. The synchrotron losses create a bump close to 
the maximum electron energies. In the Kepler remnant, the syn- 
chrotron cut-off is found to be around 0.2 keV (see Fig|2]i, while 
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Fig. 2. Energy spectrum of relativistic electrons at the shock 
front given by MHD-SDE simulations in the conditions of the 
Kepler SNR. 



in the case of G347.3-0.5 it is around 0.5 keV (see FigO. We 
verified that lowering the normalization factor qijj) of the diffu- 
sion coefficient from 16 to 3 produces a cut-off around 1 keV 
(Kepler) and 2.5 keV (G347. 3-0.5), namely that a higher cut- 
off would require a lower qip) (see Eq.[T9]). The density around 
G347.3-0.5 is badly constrained and Woo < Icm"-^ would lead 
to similar effects. We note that the above simulations maximize 
the incidence of the stochastic acceleration because we assumed 
that the resonant field dominates the total field in the downstream 
medium (see Eq.|3]l. 

To conclude it clearly appears that the downstream Alfvenic 
Mach number Vd/ VA,d cannot be much less than unity otherwise: 
(1) the X-ray filament would be too large with respect to the ob- 
served widths (see next section), (2) the X-ray cut-off frequency 
would be far larger than Ziy.cut (see Fig|2]i, and (3) the radio spec- 
trum would be harder than v (see Fig|2]i. Generally speak- 
ing, the maximum downstream resonant magnetic field cannot 
be much stronger than a few mGauss downstream of the shock 
front, otherwise regular acceleration process would be domi- 
nated by stochastic Fermi acceleration. This places an important 
constraint on the combined value of the magnetic field and the 
local ISM density as well as the respective contribution of the 
resonant and the non-resonant instability to the total magnetic 
field at the shock front. 



Fig. 3. Same plots than in Figl2]but in the G347-0.5 SNR (veloc- 
ity of the shock is 4 x \(fikmls and compression ratio rj,,^ - 4). 
The density of ISM is 1 cm"^ . The Bohm regime for the diffusion 
coefficient has been assumed with qifi - 1) = 15. The simulation 
has been performed until time t - 1600 years. 

relativistic electrons. They are displayed in Figs|4]anc|5] where 
they were obtained with parameters adapted to the dynamics of 
the Kepler and the G347-0.5 SNR, respectively. We also dis- 
played two X-ray wavebands (4-6 keV and 0.5 - 1 keV, even 
if this later wave band is usually dominated by the thermal emis- 
sion). In each case, both projected and deprojected filaments are 
reproduced. The relative normalization between X-ray and y-ray 
filaments depends mostly on the intensity of the magnetic field; 
for the same particle energy domain, it is found to scale as as 
expected. The width of the y-ray TeV rim is usually the great- 
est because an important fraction of the IC radiation is produced 
upstream. The 10-30 GeV y-rays are produced closer to the 
shock upstream than 1-3 TeV y-rays. In the downstream re- 
gion, the highest energetic electrons are confined closer to the 
shock because of their shorter radiative loss timescales. The pro- 
jected rims indicate that only a slight difference exists between 
the position of the peak of the gamma and X-ray emission. As 
the size of the y-ray rims is not much larger than the X-ray fila- 
ments, it seems impossible for any y-ray instrument to separate 
both components. This will also be the case for future instru- 
ments such as CTA unless the filaments are very lar ge (see the 
case of Vela Junior discussed in lBamba et al I (l2005ah V 



3.3.3. Comparisons between X- and y-ray filaments 

We end this section by a detailed comparison between X- and 
y-ray filaments produced by the relativistic electrons. The inclu- 
sion of neutral pion decay caused by the hadronic interaction 
with the interstellar fluid or with the shocked matter would re- 
quire a complete modelling of both the hadron spectrum and the 
ISM density profile around the SNR. This study is postponed to 
future work. 

In our calculations, the leptonic y-ray emission was inte- 
grated into two characteristic wavebands 10-30 GeV and 1- 
3 TeV using the sta ndard expression of the iso tropic inverse 
Compton emissivity jBlumenthal & Gouldlll970l) . The rims are 
produced by the scattering off the cosmic microwave photons by 



4. Diffusive shock acceleration in the case of 
downstream spatially relaxing turbulence 

We now consider a scenario where the downstream magnetic 
field fluctuations vary over a length-scale much shorter than the 
SNR shock radius /?sn- This scale noted can depend on the 
wave number k of the fluctuations. The damping of the turbu- 
lence in the downstream medium and its compression at the 
shock front can modify the particle mean residence time and the 
relativistic particle return probability to the shock. Hence, this 
magnetic relaxation is expected to modify the efficiency of the 
diffusive acceleration process itself. 

Equation ( fT2] ) shows that the particle energy spectrum at the 
shock front remains a power law, provided that the quantities 
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Fig. 4. The unprojected and projected X-ray and y-ray rims in Fig. 5. The unprojected and projected X-ray and y-ray rims in 
the conditions of the Kepler SNR (same physical conditions than the conditions of the G347-0.5 SNR (same physical condition 
Fig|2]i. For clarity, both X- and y-ray rims have been normalized than in Figl3]l. For clarity both X- and y-ray rims have been nor- 



to one. 



malised to one. 



(at a given energy E) z^ia(E) = u^idx^iAl^xiid. are large com- 
pared to unity. For Zu/d(£^) ^ 1, the particle distribution will 
be strongly softened and the acceleration timescale will shorten 
dramatically, the latter being dominated by the particles experi- 
encing the shortest residence time. A softening effect induced by 
the upstream losses is only expected at the highest energy close 
to /icR-max, namely as Zu — > 00- The diffusive length of parti- 
cles with energy lower than /scR-max is always smaller than the 
variation scale of the magnetic fluctuations (controlled by the 
highest energy), hence we have ZySE < /icR-max) 5^ 1, leading to 
a vanishing exponential factor in the above solution. In contrast, 
the softening effect downstream can be significant at energies 
much lower than ficR-max as can be highly scale (and thus 
energy) dependent. This is precisely the main topic of this sec- 
tion, namely trying to identify the parameter space that allows 
the Fermi acceleration process to be efficient in the context of a 
relaxing downstream turbulence. 

Hereafter the downstream relaxation length I a is considered 
to be energy dependent and we normalize it with respect to the 
maximum CR energy; ficR-max: 



- ^d,M X 



^CR-n 



= 4,mx|:^| . (23) 



The scale i^^A is the relaxation scale at the maximum particle 
energy 0. We reiterate recall that the relationship between en- 
ergy particle and wave vector originates from the condition that 
a given particle should resonate with a turbulence mode, krj^ ~ 1 . 
We first investigate the magnetic field profiles in the downstream 
medium resulting from various relaxation processes (Sect. 14.11 ). 
In Sect.. l4~2l the efficiency of the DSA with respect to the tur- 
bulence properties (turbulence index, relaxation index) is dis- 
cussed, by considering in particular the effect of the downstream 
magnetic field amplitude. Various numerical experiments, pre- 
sented in Section 1431 illustrate the effect of the magnetic field 
spatial variation in the particle dynamics and the associated X- 
and y-ray rims. 

4.1. Downstream magnetic field relaxation 

This work considers various turbulent magnetic field damped 
profiles: the case of an energy-dependent Heaviside profile, 
the profile produced by a non -linear Kolmogorov-type damping 
jPtuskin & Zirakashvilill2003b . and the profile produced by the 



^ Again, a correct way to handle this effect is to account properly for 
the particle back-reaction on the flow. 

^ Ail quantities with an index M are to be talcen at the maximum 
particle energy. 



Alfven or fast magnetosonic cascades (Pohl et al"2005'). We also 
briefly discuss the ca se of a turbulent dynamo action downstream 
dPelletier et alll2006h . In this section, unless specified otherwise 
5d > is implicitly assumed. 

4.1.1. Heaviside profiles 

Heaviside-type magnetic relaxation accounts for an idealized ap- 
proach to turbulence relaxation, where a given turbulence mode 
is assumed to be uniform out a distance fdik) from the shock and 
to vanish beyond that distance. This relaxation model is prob- 
ably unphysical but enables us to reproduce the basic features 
of the turbulence relaxation effects upon particle acceleration. 
Assuming this profile, we write the magnetic energy turbulence 
spectrum as (the downstream medium is defined by x > 0) 

W(k, x)d = W(x = 0+, k)n(Uk) -x) + WooTl(x - £d(k)) , (24) 

where 11 functions are Heaviside functions and x is the dis- 
tance from the shock front. The magnetic energy density far 
downstream is Wca- The normalization of the turbulent spec- 
trum W{k) = Wok^^ is related to the magnetic field at the 
shock front by means of Wo = B\x = 0+)/4;ro-A;„in, where 
^min = 2;r/l^'^^^, k = ^/Imax.d, and again cr(y6 = 1) = ln(^n,ax/^min) 
and o-(fi>l)'^ \I(J3- 1). 

The Heaviside profile, despite it crudely approximating the 
variation in the magnetic energy density downstream, permits us 
to derive a basic spatial profile of the total magnetic field given 

^-^^ ' Wik, x)dk , (25) 



At: 



which in the case of Bohm turbulence leads to (^„„„ is defined as 



< X < ^min 
tnin ^ X < ^d,M 
{d.M ^ X 



5B\x) _ 5B\0*) SBj 

An ~ An An 
5B\x) _ 6BH0^) \n{{AMlx) ^ 6BI 

An ~ An 
6BHx) ^ dBi 

An An 



Sd ln(^max/^min) 47r 

(26) 



At any given downstream location ^j.m > x > ^^in , the 
maximum non-vanishing turbulence wave number is kj^aAx) - 
kminitiMl xy^^'^ ■ Bcyond /"d.M, all turbulent modes vanish giving 
a total magnetic field B^o close to the ISM magnetic field value. 
The spatial variation in the magnetic field for any other diffu- 
sion regime is more complex, as it scales as 1 - (^d.ivi/j^)*'^^^'*'' 



Marcowith & Casse: Postshock turbulence in young supernova remnants 



11 



for X > ^inin. The total magnetic field is required to calculate the 
synchrotron losses properly, in addition to the normalisation en- 
tering the particle Larmor radius and the local Alfven velocity. 
Once the total magnetic field is known, we can calculate, for ev- 
ery relativistic particle of energy E, is the fraction of the total 
magnetic field that can resonate with this particle, namely by in- 
tegrating all turbulence modes verifying 1 /?/,(£■) < k < femaxC-^)- 
This is achieved by computing the function b defined in Eq.(l7]i. 
If magnetic turbulence relaxation follows a Heaviside prescrip- 
tion then one obtains 



b(Q<x< i„ 



^ 4oh I Mm 



(27) 



^coh 



hXE) 

^coh 



CR-maxJ 



''coh 



Once both the total magnetic field and function b are known, it 
is easy to compute in our simulations both the spatial and energy 
diffusion coefficients of every test particle, which are mandatory 
to determine the particle motion and stationary particle distribu- 
tion solutions in Eqs. ( fTTT i and ( IB. 21 ). The procedure is repeated 
in the same way for any magnetic profile. 

4.1 .2. Non-linear Kolmogorov damping 

In models of incompressible MHD turbulence described by 
the Kolmogorov energy cascade towards large wave numbers, 
the non-linear damping kernel scales as k^^^W{k)^^^. Following 
iPtuskin & Zirakashvilil (l2003h . this kernel can be simplified 
while still respecting the spatial relaxation profile. We have 



rNL(/fc, x) ^ To X k^'^W(k, x) 



1/2 



(28) 



where To - 5 x 10"^ x V^^/iBj/Any'^. Here we consider the 
cascade to be initiated behind the shock and use the local total 
magnetic field and Alfven velocity. 

In the shock rest-frame, the turbulence relaxation downstream 
(for X > 0) is described by a stationary equation 

Vsh dWjk, x) 

— x = -2rNL(^, x)W(k, x) , (29) 

rtot ox 

and a boundary solution W(k, x = 0^) = Wo x k'^. The solution 
ofEq.(|29llis 



W(k, x) 



W(k, X = 0+) 

Xq 



(30) 



An estimate of the scale xq is fsee lPohl et all (120051) -) 

1/2-1/2-1/2 

xo-K - [300 X i„„,d] X 



'"tot 



(31) 



We used & - cr/ 16 and the shock velocity ysh,4 is expressed in 
units of 10"* km/s. The downstream maximum turbulence scale 
^mdx,d can be connected to the maximum Larmor radius of CRs 
upstream by means of Eq.(fT4l). Reduced rigidity at maximal en- 
ergy EcR-mix is such that pLu - Pm as the diffusion coefficient 
rapidly increases as E^ beyond /icR-max- Both conditions set the 
maximum upstream turbulence scale /tmax.u and the maximum 
CR energy ^cR-max- We find that .Imax.d - 5.2 rL-max,u/''subPM, 
where pM - Pm/0.3. 

The relaxation scale is fd(E) - {d,M x {E I Ecu-maxP^^^^^ ■ 
The factor {^jsa is defined as the length over which turbulence 



level has decreased by 1/e compared to its value at the shock 
front, i.e., ^d,M = ( V^- 1) ^0- The spatial dependence of the total 
magnetic field and function b were calculated using Eqs. ( 1251 ) 
and (IT). These expressions, which are quite lengthy especially 
for the b function in Eq.(|6]), were implemented into the code but 
are not explicitly given here. 



4.1.3. Exponential profiles 

When turbulence damping rate does not depend on space but 
remains dependent on wave number (F = T{k)), the relaxation 
of the downstream magnetic field follows an exponential cut-off 
on a scale length ti{k) - rioJ'{k)IVsh. The turbulent magnetic 
energy spectrum is then 



W(k, x) = W{k, 0+) X exp 



(32) 



The Alfven and Magnetosonic waves cascades considered by 
iPohl et ai (l2005h follow this scaling, the corresponding damp- 
ing rates and expression for xo can easily be obtained from their 
Eqs. (8) and (11) respectively. Considering the Alfvenic cascade, 
we obtain 



X<3-A 



T/ 1/2-1/2 



tot 



d,-4 



(33) 



The coherence scale of the downstream turbulence is /"coh = 
'imax.d Pm/27!'. The fast magnetosonic cascade leads to a similar 
expression except that the wave phase velocity can be approxi- 
mated as VpM.d - (^Ad ''sd-''^^' "-s-d being the sound velocity 
behind the shock front. We note that the above expression for the 
Alfven cascade results from the combinaition of the critical bal- 
ance and the anisotropy obtained i n the Goldreich-Sridhar phe - 
nomenology of strong turbulence ( Goldre ich & Sridhar 1 11991 . 
Again, the expressions of the total magnetic field and resonant 
field are rather lengthy and are not shown here. We note that 
Eqs. ( ISTb and (|33T ) show that the Kolmogorov damping leads to 
a slower cascade and thus to longer relaxation scales than expo- 
nential damping. 

4.1 .4. Turbulent dynamo downstream 



iPelletier et a] (l2006h (see also IZirakashvih & PfiiskinI (l2008h ) 
discussed the action of a turbulent dynamo in the downstream 
medium that would lead to additional amplification of the mag- 
netic field. The magnetic field is expected to saturate at values 
close to equipartition with the dynamic gas pressure. The dy- 
namo action is driven by the non-vanishing helicity of the non- 
resonant turbulent modes. 

The coiTesponding scale of magnetic field variation is given 
by the ratio of the magnetic turbulent diffusivity Vt to the dy- 
namo amplifi cation coefficient ap . The two coefficients can be 
expressed as (IPelletier et alll2006l) 

2c / V' \^ 

ttD - ^ X ( -^j X In (rL (ficR-Max) //"L (^CR-min)) , (34) 



and 



2cA„ 



37r2 



Vsh 



(35) 



where Zscr min represents the lowest resonant energy. The ampli- 
fication scale is then (am^\ ~ A^axA/i^fp)- Turbulence modes of 
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wavelength longer than {^mpi grow and saturate close to equipar- 
tition. Other turbulence modes are expected to dampen rapidly 
(over a few plasma skin depths) because the non-resonant waves 
are not normal modes of the plasma, as already stated in sect. 
2:21 



4.2. Particle acceleration in a relaxed-compressed 
turbulence 

In the next few paragraphs, we present some useful analytical 
estimations for the analysis of the numerical simulations pre- 
sented in sect. 14. 31 These calculations used the Heaviside related 
profiles derived in sect. K.l.D We note that the following charac- 
teristic timescales are strictly valid in the framework of infinitely 
extended diffusive zones but are used to discuss the effect of a 
spatially limited diffusive zones. However, we see in sect. 14.31 
that these approximations lead to correct energy spectrum fea- 
tures, except at the highest energies. 

4.2.1. General statements about turbulence parameters 



iPohl et all (l20Q5h discussed various possible downstream relax- 
ation processes. First, the non-linear Kolmogorov damping pro- 
duces a relaxation length idik) oc k^-^)l-. Each turbulence mode 
k being in resonance with relativistic particle whose Larmor ra- 
dius verifies k?]^ > 1, we obtain 6^ - (3 -y6)/2 > (between 
1 and 1/2 for 1 < ^ < 2). The two other processes considered 
bv lPohlet all (12005) scale as k namely 6^-1 /2. A variation 
range of 6a between 1/2 and 1 is then clearly identified. We ex- 
tend it to encompass the regime 6^ - Q, a. limiting case where 
relaxation lengths are spatially independent. 

What if 6i were negative ? A strict lower limit to 6^ is given 
by the condition ^d(£cR-min) ^ ^sh- A non-relativistic minimum 
resonant energy -EcRmin - 0.1 x ( V2 - l)mpC^ seems accept- 
able so that 5d > 5d,lim = ln(/?sh/^d,M)/ln(£'cR-min/£'cR-max)- 

The lower limit 5d,iim has typical values of between -0.3 and - 
0.2 when identifying ^d.M with the size of the X-ray filament. 
Relaxation regimes with Sd < Q do not necessary correspond to 
any known damping process but have some interesting proper- 
ties, in particular concerning the radio filaments. 

4.2.2. The dominant loss mechanism 

Comparing typical energy loss timescales is a useful tool for 
determining whether or not diffusive particle losses can affect 
the energy spectrum of relativistic particles. Assuming that tur- 
bulence relaxation follows a Heaviside prescription, we can ex- 
press these timescales by assuming a constant downstream mag- 
netic field on the relaxation length relative to a particle of 
energy E. 

Four timescales are relevant to set the maximum particle en- 
ergy in a relaxed and compressed turbulence: 



1. The acceleration timescale is given by 



Pd(E) \' 



^m3.x d — 2 

facc(£') - [7 yrs] X pm X y(r)K(y6, r) x — x I , 

(36) 

where K(j8, r) - q(J3) x (H(J3,r)/r + I) and the maximum 
wavelength of the downstream turbulence is expressed in 
units of 10"^ pc. 
2. The advection timescale, i.e., the time required for a particle 
to travel over a distance while being advected with the 



downstream flow, is given by 



taiy(E) 



Vd 



[4yrs]xrV-,yd-2(E). (37) 



3. The diffusive timescale, i.e., the time required for a particle 
to travel over a distance in a diffusive motion. Q 



tdiff(E) = -TTTTT^ - [0.3 yrs] X x 



2 /pd(E)\^-^ 

max,d-2 \ PM j 



(38) 



6Dd(£) 

4. The synchrotron loss timescale 

fsy„(£) - [1.25 X 10^ yrs] x E^iv x B,!^ x f,j,„c , (39) 



where the parameter fsync is represented by {Hifi, r) + 
r)/(H(J3, r)/ + r). This expression takes into account the 
mean residence time both in the upstream and downstream 
medium. 

The maximum electron energy is given by the equality 
facc(£'e-max) = floss (-Ee-max), whcrc floss is the shortcst of the syn- 
chrotron, advective, and diffusive timescales. When X-ray fila- 
ments are controlled by the radiative losses, we have fioss = fsyn- 
In the case of escape losses being the most significant losses 
for filaments, we then have fioss = min(fdiff, fadv)- It can be 
seen that particles of energy close to Ee-max, diffusive losses 
are always dominant compared to the advection losses, hence 
fioss(£^e-max) = fdiff (^e-max)- We notc that the downstrcam resi- 
dence time fi-es.d - ( Vd/c)facc (during one Fermi cycle) should not 
be compared with diffusive or advective timescales because only 
particles returning to the shock experience a full Fermi cycle. 
Performing this comparison would lead to a maximum particle 
energy /Se-max much higher than values obtained in the context 
of our numerical simulations. 



4.2.3. Conditions for an efficient particle acceleration 

For relaxation-dominated filaments, the ratios of the accelera- 
tion timescale (Eq.[36]l to either the diffusive (Eq.[38ll and to the 
advective (Eq.l37Ti timescales vary as £'2(2-/3-5d) ^jj^j j7(2-fi-Si)^ j.g_ 
spectively. Two different regimes are now discussed. 



2-Sd-/3>0: Once E < E, 
< fdiff and t, 



e-m-dx, the various timescales order 
the acceleration process can oc- 



as face < fdiff and face < fadv 

cur without noticeable losses and thus a particle energy spec- 
trum behaves as a power law. We note that for energy lower than 
Eadvj advection losses become dominant compared to the dif- 
fusive losses. Formally, we derive this energy limit by setting 

racc(£'e-inax) = fdiff (^e-max), which Icads tO 



adv 



6 



l/2(2-5j-/5) 



(40) 



where g{r) = 3/(r - 1) x {H(J3, r)/r + 1). We hereafter note fliat 
sir) = g(r)/g(4). 

2 - Sd - p < 0: In this case, the ratio of the diffusive to advec- 
tive timescales is always lower than unity, i.e., diffusive losses 
dominate at all energies. Once E < iie-max, downstream escapes 
limit the shock acceleration process considerably as the accel- 
eration time becomes longer than fdiff as energy decreases. The 



' The factor 6 in the denominator of Ea.ll38t. which appears to be the 
random walk along the radius of a sphere is composed of 3 independent 
random walks along each cartesian coordinates 
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SNR 


^d-diff 




^d-(iml ^d-diff 


^d-Uml ^d-diff 


q(l) = 1 


q(i) = 1 


q(\) = 1 


q{i) = 16 


Cas A 


311 


394 


2.3 


0.4 


Kepler 


220 


293 


3 


0.5 


Tycho 


210 


333 


5.1 


0.8 


SN 1006 


174 


189 


1.2 


0.2 


G347.3-0.5 


164 


183 


0.95 


0.15 




5d =0 


6i = 1/2 


V(5d 


V<5d 



Table 3. Table presenting analytical estimates of the downstream magnetic field valu e in the context of diffusive-loss-dominated 
SNRs rims. The SNR rim observed parameters are the same as in iParizot et all (l2006l) and the shock compression ratios are rtot = 

'"sub = 4. 



same conclusion can be reached from a close examination of the 
particle distribution given in Eq. (fT2l i. The term Zd = u^td/Dd is 
proportional to (fdiff/4cc)'^^, and 2- jS- 5a < leads to Zd tend- 
ing toward zero. The particle energy spectrum then steepens at 
low energy, which is obviously in complete disagreement with 
the Fermi acceleration scenario. 

Hence, efficient Fermi acceleration is only possible if 2 - 
6a —p > 0. For instance, an energy independent relaxation length 
5i - (as well as 6^ < 0) verifies this criterion for all diffusion 
regimes. In the case of a Kolmogorov type non-linear turbulence 
damping, the supplementary relation 5i - {?> — y6)/2 imposes 
/3 < 1, which means that only the Bohm regime can fulfil the 
previous condition (we see in sect. l4.3l that particle acceleration 
is inefficient in that case). In the context of Alfven and magne- 
tosonic cascades, Kolmogorov turbulence regime (j3 = 5/3) is 
the sole regime faiUng to verify the previous condition. 

4.2.4. Magnetic field limits in a relaxed-compressed 
turbulence 

In the context of X-ray filaments controlled by the down- 
stream turbulence damping, we can link the size of the fil- 
ament, noted ARx, to the maximal relaxation length ^d,M as 
^d,M = A^x(£'e-max/£'e-obs)'''' = C{6d)tJix- El The energy Ee-obs 
is the energy of particles emitting in the 4-6 keV band and this 
value depends on the local value of the total magnetic field. 

A downstream magnetic field estimation Bd.diff can be ob- 
tained from the dynamics of the electrons by requiring that 
4cc(£'e-max) - fdiff(£'e-max) using the prcvious relation between 
and hRx- In the context of the Bohm diffusion, one obtains 



fid,-4,diff ^ 3.7 X ^oe = 1) 



2/3 



(E 



1/2 

y-cut,keV 



y(r) 



1/2 \ 



A/?x,-2C((5d)Vd,- 



2/3 



(41) 



where Vd.a = VijlQ^ km/s and again g(r) - g{r)/g(r - 4). If 
/? > 1, the derivation of the magnetic field amplitude is more 
cumbersome. 

The determination of the SNR X-ray filaments are domi- 
nated wether either by the relaxation of the downstream mag- 
netic turbulence or synchrotron losses is provided by the con- 
dition fdiif(£'max) = fsyn(£'max) = facc(£'max)- The Corresponding 



* The dependence of on the wavelength /I is a priori valid only up 
to /Imax - ^L(£cR-max) and. We should strictly not expect the scaling of 
^d to extend beyond A^^^. Above A^.^^, the diffusion coefficient increases 
as R^^ and particle acceleration continues to proceed beyond £cR-max but 
the number of particle accelerated and the turbulence energy density 
both rapidly drop. For this reason, we consider 6d to be controlled by 
the kernel of the damping rate above -BcR-max as, e.g., in the case of the 
Kolmogorov damping 6i = 3/2 in this energy regime. 



limit value of the magnetic field is (again in case of Bohm diffu- 
sion) 



fid,-4,lim - 8.9 X 



* sh,4 7-cut,ke V „ 



,2/3 



(42) 



To ensure that SNR X-ray filaments dominated by the relaxation 
of the magnetic field, it is compulsory to have Bd.diif < fid.iim- 
The factor qifi) was isolated in Eq. ( 1411 1 to show that no solution 
is then possible if qifi = 1 ) s> 1 . In other words, a diffusion coef- 
ficient close to the Bohm value is required to allow the relaxation 
of the turbulence to control the size of the filaments. We have 
also to keep in mind that the downstream magnetic field am- 
plitude has to be consistent with the aforementioned assumption 
that an amplification upstream has occurred, namely » Bism- 
Concerning cosmic rays of energy E ~ EcR-mux, the down- 
stream diffusive losses will dominate if particles cannot es- 
cape from the upstream region into the ISM. This imposes a 
constraint on the magnetic field amplitude at the shock ob- 
tained from Eq. (l22l i. Indeed, upstream escape losses are dom- 
inant if facc(£'cR-max) < fditf (^CRmax), whcn wc assumc that 

^d(£'cR-max) - A/?x(£'cR-max/£'e-obs)'''' ■ 

- In the case 6d - 0, downstream diffusive spatial escape 
downstream always controls the maximum CR energy. 

- For Sd + 0, the previous condition leads to an upper limit to 
the downstream magnetic field, noted Bd.esc- Hence if > 
Bd.esc the CR maximum energy will be fixed by the upstream 
escape losses and conversely if B^ < Bd.esc the CR maximum 
energy will be set by the downstream escape losses. 

The downstream magnetic field then has to fulfil Bd.diff ^ 
min(Bd,esc, Bdjim) to ensure that the downstream turbulence re- 
laxation be the controlling process of the energy cut-off of rel- 
ativistic particles. Applying the previous conditions to our SNR 
sample, we always find that Bd.iim < Bd.esc- This means that an 
intermediary regime may exist where electrons lose their energy 
through radiative losses while cosmic rays cut-off is set by down- 
stream diffusive losses. Of course, if the magnetic amplification 
process is efficient enough to generate higher turbulent magnetic 
field amplitude then upstream losses will take over. 

Table (|3]l displays the values of Bd.iim and Bd.diff related to 
our SNR sample. The Kolmogorov regime was discarded as it 
does not produce any efficient acceleration as we see in sect l4.3l 
We show that for SN 1006 and G347. 3-0.5, not much room is ef- 
fectively left for the case of magnetic relaxation-controlled fila- 
ments. This result seems rather robust as a variation of the shock 
velocity by a factor of 40%, or a variation in the synchrotron 
cut-off by a factor of 2 does not lead to any variation in the mag- 
netic field greater than 25%. However, a variation in the filament 
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Electron spectrum at front sliock Acceleration tinicsealc 




Distdiioe from tlie shock (pc) Distance from tlie sliool; (po) 

Fig. 6. Energy spectrum of relativistic electrons at the shock 
front given by MHD-SDE simulations in the conditions of the 
Kepler SNR (see Figl2]for details). 

width by a factor of 2 would imply a variation in the magnetic 
field by a factor of 60%, which may slightly modify the previous 
conclusion. Quite generally, the maximum magnetic field ampli- 
tude is found to be in the range ~ 200 - 300 /iGauss. 
To summarise we find that if downstream magnetic relaxation 
controls the features of the SNRs X-ray filaments, a Bohm-like 
diffusion regime is likely to occur and the particle diffusion co- 
efficient normalization factor qip - I) has to be quite close to 
unity, i.e., the diffusion regime has to be close to a genuine Bohm 
diffusion regime. In this context, we show that only a fraction 
of our SNR sample to meet these conditions, namely the young 
ones. Using the various observational constraints related to the 
older SNRs (SN1006 and G347.3-0.5), we have shown that the 
X-ray filaments existing in these objects are likely to be ruled by 
radiative losses associated with synchrotron emission. 



ation, the stochastic reacceleration, and the radiative losses for 
the electrons. In the following paragraphs, we discuss the phys- 
ical agreement between assuming magnetic field relaxation to 
control the X-ray filaments and the actual results coming from 
the computation of relativistic electrons acceleration. 

4.3.1. Downstream magnetic Kolmogorov damping 

When non-linear Kolmogorov damping occurs in the down- 
stream medium of the shock, we have seen in the previous sec- 
tions that two conditions have to be fulfilled to reproduce both 
the appropriate energy cut-off and the correct size of the ob- 
served X-ray filament. These two conditions can be expressed 
as: having the correct downstream magnetic field given by 
Eq.(l4TI) (to ensure that the electron energy cut-oflF is consistent 
with the observations) and having the typical magnetic relax- 
ation length xo-K (see Eq.lSTTi that is similar in the size of the X- 
ray filament. In the non-linear Kolmogorov regime, the only dif- 
fusion regime able to provide to an efficient particle acceleration 
is the Bohm diffusion regime, where the relaxation energy index 
Sd - ^- Inserting, for the Kepler SNR, this value into Eq. (14111 
leads to a downstream magnetic field of ^ 390//G and a re- 
laxation of xo-K - 0.39 pc. The relaxation size is clearly too 
large to provide an X-ray filament, whose thickness is inferred 
to be of the order of 10"^ pc from X-ray observations. Applying 
the same reasoning to the other SNRs leads to a similar conclu- 
sion: having both the appropriate electron energy cut-off and X- 
ray filament size is incompatible with a non-linear Kolmogorov 
occurring in the downstream medium of the SNR shock. The 
only way to overcome this conclusion would be to have the fac- 
tor cr - ln(femax/^min) to be much smaller than expected (see 
EqETJ. Anyway, having cr so low would mean that the range 
of particle energy able to resonate with turbulence mode would 
be so narrow that it would not be able to provide any signifi- 
cant acceleration. This explains why our result differs from the 
conclusion drawn by Pohl et al (2005). It seems then that it is 
very unlikely that non-linear Kolmogorov damping, which is a 
slower process than Alfven/ fast magnetosonic cascade, occurs 
in the downstream medium of SNR shocks. 



4.2.5. Radio filaments 

The energy of the radio electrons is typically four order of mag- 
nitude below that of the X-ray emitting electrons, i.e., 

£eobs,R-[1.5GeV] 

where iSy-obs-GHz is the energy of the radio electrons emitting 
in the GHz band. Using both Eqs. (|39] | and dJTl l. one can eas- 
ily check that the synchrotron loss timescale at iieobs.R is always 
longer than the advective loss timescale, unless 6d is lower than 
typical values of the order of -0.5, a value always lower than 
^d.iim- If <5d-iini < ^ 0, the Small turbulence scales relax on dis- 
tances longer than ARx- This very particular case would produce 
radio filaments larger than the size of X-ray filaments inferred 
from the Chandra observations. In contrast, the regime 6^ > Q 
would allow the largest fluctuating scales controlling the size of 
the radio filaments. In this cas e, the radio filaments a re expected 
to be of the order of ARx (see lCassam-Chenai et al.l (i2007i) ). 

4.3. Numerical simulations 

We performed MHD-SDE simulations by taking into account all 
previous settings, namely the downstream magnetic field relax- 



4.3.2. Alfvenic-fast magnetosonic mode damping 

In the context of Alfvenic-fast magnetosonic turbulence relax- 
ation, the typical relaxation length xq-a is shorter than xq-k- 
Compiling the aforementioned necessary conditions to repro- 
duce accurately an X-ray filament in the SNR environment, we 
obtain a typical xq-a of the order of 10"^ pc when using mag- 
netic field values provided by Table(|2]l. This means that the 
Alfvenic-fast magnetosonic modes damping is a plausible candi- 
date to explain the presence of SNRs X-ray filaments. To sustain 
this conclusion, we performed, in the context of the Kepler SNR, 
MHD-SDE simulations designed to reproduce the dynamics of 
relativistic electrons and the associated X-ray and y-ray emission 
maps. In figures |6] and |7] we display the particle distribution at 
the shock front and the X- and y-ray filaments respectively. All 
simulations were performed in the Bohm regime. In that case, 
2-Sd - p - 1/2 > 0. In each cases the magnetic field is damped 
in the downstream medium following an exponential relaxation 
as in Alfvenic-fast magnetsonic modes damping. Bohm regime 
in downstream region has been assumed. The dashed-line shows 
the stationary solution found in Marcowith et al (2006), which 
includes particle reacceleration in the Fermi cycle. In the up- 
per right panel the acceleration (only the regular Fermi accelera- 



Marcowith & Casse: Postshock turbulence in young supernova remnants 



15 



Deprojected 7 and X — ray r 



B„=293 


7-rcy 




(1-3 TeV) _ 


N,s.=0.7cm"' 






-■- 




1 (10-30 GeV)" 




i (4-5 keV) 








t\ (0.5-1 keV)- 







Projected 7 anc| X — rqy 1 



V,„,=5400 km/s 


7-roy ' 




Ts,p=400 yrs 


-.- 7- ray 




(10-30 GeV) 




X-ray 




(4-6 keV) 




X-ray 




(0.5-1 keV) " 




1 

\ 



Electron spectrum at front sliock 



Acceleration timcscalc 



-0.02 -0.01 0.00 0.01 0.02 
Distance from the Shock (pc) 



-0.2 -0.1 0.0 0.1 0.2 

Distance from fhe Shock (pc) 



Fig. 7. The unprqjected and projected X-ray and y-ray rims in 
the conditions of the Kepler SNR in the case of an exponential 
relaxation profile. For clarity, both X- and y-ray rims have been 
normalised to one. 



tion), and the diffusive and downstream residence timescales are 
displayed using solid and dashed lines. Diamonds represent for 
our numerical calculation of the acceleration timescale, which 
is in agreement with the theoretical estimation. We also display 
in the two lower panels the spatial dependence of the diffusion 
coefficient at the maximum electron energy (lower left) and the 
magnetic profile in the downstream medium at f = 400 years 
(lower right). 

Several obvious differences appear in both figures |6] and [T] 
with respect to the simple advection case presented in figures |2] 
and [3] First, as stated in sect. 14.21 the normalization of the dif- 
fusion coefficient qiJS) has to be close to one. Even in this case, 
the maximum particle energy is limited to values close to ten 
TeV (for parameters associated with the Kepler SNR). One of 
the necessary conditions to fit the observed size of the X-ray 
rim, namely xq-a ~ ARx, produces an increase in the diffusion 
coefficient by a factor of a few tens above the typical diffusion 
length, and consequently low maximal energies for both elec- 
trons and cosmic rays. The X- and y-ray filaments also exhibit 
some different features in the case of an Alfvenic-like relaxed 
turbulence. The low energy particles producing the synchrotron 
photons in the interval 0.5-1 keV and the y-ray photons in the 10- 
30 GeV band, respectively, do extend to shorter distances behind 
the shock (electrons having energy ~ 1 TeV). This can be under- 
stood by the effect of the resonant component of the magnetic 
field b in Eq.®. At a given downstream location, particles with 
energies E < fimax do interact with a lower number of modes 
than in the advected case. This effect is caused by high wave 
number modes relaxing over shorter distances than lower wave 
number modes within the same turbulence spectrum. Compared 
to the advected case, more low energy particles experiencing dif- 
fusive losses are lost than at highest energies (which are also 
subject to diffusive losses). Particles of energy around a few tens 
to hundreds of GeV are then confined to closer to the shock and 
do not experience strong magnetic field variation: the standard 
shock solution is then recovered in this domain. We verified that 
the shock synchrotron spectrum cuts off at an energy close to 
one keV. 

We also tested the solution in the case /3 - 2, i.e., 2-6d -/3 - 
- 1 /2 < 0. No significant particle acceleration has been found as 
diffusive losses dominate at low energy (see FigO. The numer- 
ical acceleration timescale is also found to be shorter than the 
theoretical estimation which is consistent with particles only re- 
turning quickly to the upstream medium after entering the down- 
stream region, are able to avoid massive diffusive losses. These 
simulations confirm the conclusions drawn in sect. 14.2.31 
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Fig. 8. Same case as treated in Fig|6]but with P - 2. Here, mas- 
sive diffusive losses are occurring since 2 - 6d - /3 < and thus 
no significant acceleration is observed. 

4.3.3. Solutions for turbulent dynamo amplification 

The coherence length of the downstream turbulence entering 
the evaluation of ^ampi in sect. 14.1.41 cannot be longer than the 
X-ray filament width, otherwise the condition about the maxi- 
mum CR diffusion coefficient upstream given by Eq. (l22l i would 
not be satisfied. This means that if a magnetic dynamo oper- 
ates downstream, then the growth scale length is < ARx- The 
growing modes are restricted mostly to large scales, i.e., to wave 
numbers close to kj^in- They are considered for the particles to 
contribute to the mean magnetic field. The rapid increase in the 
magnetic field downstream to values close to equipartition pro- 
duces enhanced radiative losses and thus much thinner filaments. 
We checked the effect by performing simulations in which we 
added a mean magnetic field downstream of values close to a 
few mGauss. 



5. Discussion and summary 

Young SNRs are strong particle accelerators, as illustrated by the 
presence of thin X-ray filaments. In these astrophysical objects, 
the X-ray emission is produced by synchrotron radiation, involv- 
ing particles whose maximal energy is higher than tens of TeV 
and mag netic field strengths behind the shock of a few hundred 
/iGauss (Parizot et al 2006) . This work has extended the study 
undertaken bv iParizot et all (l2006h of the physical properties of 
both the turbulence and transport coefficients in the same sam- 
ple of five young SNR. We have included the turbulence com- 
pression at the shock front, the possibility of particle reaccel- 
eration in the downstream region of the shoc k, and the relax - 
ation of the magnetic fluctuations downstream (iPohl et alll2005h . 
We have also described the generation of magnetic fluctuations 
in the shock^precursor for the two regimes of the streaming in- 
stability (iPelletier et alll200 6 |). This work has been developed in 
the same framework as iLa gage & Cesarskvl (Il983h but adapted 
to the case of amplified magnetic fields around SNR, although 
the maximum CR energy has not been fully investigated here. 
We have developed a numerical scheme based on the coupling 
between the equations of magnetohydrodynamics and a kinetic 
scheme handling the calculation of the electron particle distribu- 
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tion function. The scheme involves a set o f stochastic differen- 
tial equations (SDE) described elsewhere (ICasse & Marcowi"5il 
[2003, 2005). The SDEs have been adapted to account for the 
discontinuity in the diffusio n coefficients properly using a skew 
Brownian motion (see also IZhangI (l2000h ). The following con- 
clusions can be made: 



- The normalization (i.e., factor qifi)) of the spatial dif- 
fusion coefficient should remain close to unity to avoid 
massive particle diffusive losses, leading to a drop of the 
Fermi acceleration efficiency. A genuine Bohm diffusion 
regime is then required if magnetic turbulence relaxation 
is to occur in the downstream region of the shock. 



1 . The compression of turbulent scales at the shock front does 
not deeply modify the efficiency of shock acceleration. The 
conclusions addressed by Parizot et al (2006) are found to 
be robust in the case of a downstream, advected, magnetic 
field, young SNRs exhibiting X-ray filaments do accelerate 
particles to at most PeV energies. 

2. For the various regimes of streaming instability occurring 
in the shock precursor, the SNRs contained in our sample 
are expected to generate magnetic fields up to a few hun- 
dred /iGauss. For shock velocities of a few hundred thou- 
sand km/s, the level of fluctuations tends to be shared by the 
non-resonant and the resonant regimes. The resonant modes 
may contribute to some particle reacceleration downstream. 
However, the amount of reacceleration cannot be too large, 
otherwise the shock particle spectrum would be harder and 
the X-ray filament width would be larger than observed. This 
provides an observational constraint of the number of reso- 
nant modes present downstream of the shock front. The fate 
of non-resonant modes generated upstream still requires con- 
sideration. 

3. We have presented calculations of the projected and depro- 
jected X- and y-ray filaments, each one in two specific wave- 
bands. For the separation between the X and y-ray peak 
emission is found to be far below any y-ray mission resolu- 
tion capabilities when observing young SNR, some detailed 
observations could be undertaken for more extended objects 
such as Vela Junior. 

4. For relaxed turbulence occurring in the downstream region, 
our conclusions are the following: 

- When the magnetic relaxation scale varies as tiik) oc 
k^^'', a magnetic turbulence (whose power-law index is 
P) is able to provide suitable conditions giving rise to an 
efficient particle acceleration \f2 - 6^ - fi > Q. 

- We have tested several relaxation processes obtaining 
various values of 6^. When Kolmogorov damping occurs 
in a Bohm diffusion regime, it appears unlikely to pro- 
duce strong acceleration in the framework of relaxation 
limited filaments when accounting for the complete dy- 
namics of the turbulent spectrum. On the other hand, the 
Alfven and fast magneto-sonic cascades provide suitable 
conditions giving birth to particle acceleration while be- 
ing able to match all observational features of X-ray fila- 
ments. In this context, we have found that the maximum 
energy particle (both for electrons and cosmic rays) can- 
not be much higher than a few tens of TeV. 

- The magnetic field strengths downstream of the shock 
cannot be much higher than 200 - 300/v Gauss, otherwise 
radiative losses would control the X-ray filament width. 

- For the supernova remnants SN1006 and RXJ 1713- 
3946.5, none of the various turbulence relaxation pro- 
cesses considered in the present paper have been able 
to provide efficient particle acceleration and match the 
corresponding observational features. It seems that only 
the youngest SNRs (Tsnr < 500 yr) of our sample may 
exhibit X-ray filaments controlled by downstream turbu- 
lence relaxation. 
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Appendix A: Magnetic field profile produced by the 
resonant instability 

The amplification factor related to the resonant instability de- 
pends on t he amplification facto r produced by the non-resonant 
instability (IPelletier et all (l2006h . Eq.34) and is given by 

A\{x) - Q-res X Anr(.x) X I d\n{k) (exp(-fl(jc)^^"^) - 1/ej , 

(A.l) 

where aj-es - tt/i/ixMacxj^cr > 1 and is the maximum resonant 
wave length at a distance x, k - kt^oh varies between feinin(= 
l/'-L(£'cR-max))4oh - 1 and ^,(x)4oh > 1 0. We have: 



a(x) = — X (Vsh/c) X (x/4oh) X /7tot(x) < 1 
/3(f> 



The exact integration of Eq. dA.ll i involves a difference between 
two exponential integral: Ei{-a(x)kt) - Ei(-fl(x)). The second 
term dominates when k, > I, and we obtain 

Ar(x) cc [Anr(x) X {-Ei(-a(x))/(2-p) - \n(k4x))/ sxp(l))Y'^ . 

(A.2) 

The above equation is implicit because the total magnetic field 
is hidden in A;, and ritot- 

At distances x <s: ^diff(£^CR-raax) where a(x) <k 1, we approx- 
imate -Eii-aix)) ^ - ln(a(x)) - C, C - 0.5772 is the Euler 
constant. At a first approximation, within the precursor Ar(x) 

1 /2 

scales as A(x)^^. 



Appendix B: Derivation of the shock particle 
distribution function 

The steady-state general ID Fokker- Planck equation is given by 



df d I df 
u— - — \D— 

ox ox \ ox 



+ (Md - llu)6{x) 



df 
d In 



(B.l) 



Where the upstream medium is defined by -iuip) < x < 
and the downstream medium by < x < £d{p)- The shock 
front is at X = 0. In this equation, we have neglected the syn- 
chrotron/turbulence generation losses since we focus on the par- 
ticle diffusive losses. The presence of finite extensions in both 
the upstream and downstream media imposed by boundary con- 
ditions for / as f{-(u7 p) = = fifdi p)- To determine the spatial 
behaviour of the / function, we integrate Eg. (IB. lb from the left 
boundary to x in the upstream medium and from x to the right 



' As discussed in sect. 12.1.11 we assume the same coherence length 
over the whole precursor. 
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boundary in the downstream medium, and we obtain 
J^^ expij^^ 6u(x" , p)dx")dx' 



Mx,p) = fsip) 



fd(x,p) = fs(p) 



exp(£^^^ Ouix' , p)dx')dx 
J^" exp(- J^' d(x",p)dx")dx' 
J^^' exp(- J^' ddix', p)dx')dx 



(B.2) 



where /s is the distribution function evaluated at the shock front 
and the functions Quid are the inverse of the effective diffusive 
lengths and defined to be 



9u/d(Jc, p) 



dx 



(B.3) 



u/d 



The energy flux carried by the relativistic particle has to be con- 
served throughout the shock front, namely for i/ — > 



dx d In p^ 



= 0. 



(B.4) 



The spatial derivatives of / are evaluated using Ea.( IB.2b . which 
produces a differential equation for /s: 

d\n fs (p) 3 I ^u^O' ^'^P^Xl ^1*^-^'' Z')'^-^') 

dlnp (u^-Ui) I £° exp(£^ 0u 

^ Z)d(0,p)exp(- J^gd(x^pW) ] 
j^^" exp(- f^' Odix' , p)dx')dx j 



(B.5) 



Appendix C: Particle acceleration and multi-scale 
simulations 

This section presents the numerical framework used to simulate 
both the supernova thermal plasma evolutio n and the relativistic 
charge d particles transport. As detailed in ICasse & MarcowiflJ 
(l2003h and |Casse & Marcowith (20051), the background fluid 
and large-scale magnetic field are calculated using the m agne- 
tohydr odynamics code VAC for Versatile Advection Code (iTothI 
(119961) V The simulations are performed using a ID spherical 
symmetry, where the evolution of the supra-thermal electrons 
and nuclei are calculated using t he stoc hastic differential 
equations (SDE) formaUsm (iKriiUs & Achterbergi Il994t) . The 
numerical description of supra-thermal particle transport is 
crucially dependent on the ability of the MHD code VAC 
to capture the shock structure. To obtain the sharpest shock 
front possible, we used the TVD-MUSCL scheme coupled 
with a Roe-type approximate Riemann solver dToth & Odstrcill 



Section iCTI brieflv reports on the MHD-SDE schemes used 
to model a ID spherical SN remnant expansion. In particular, 
sects IC.2I and IC.2.11 discuss at length the stochastic differential 
Euler schemes with spatially dependent diffusion coefficients 
and their application to the diffusive shock acceleration prob- 
lem. Section lC.2.2l describes the shock capturing procedure that 
efficiently couple the MHD and SDE schemes. 



C.1. Supernova remnant modelling 

The time evolution of the thermal magnetised plasma is fully 
controlled by the MHD equations providing mass, momentum. 



and energy conservation as well as electromagnetic field induc- 
tion, namely 



§ + V ■ (pV) = . 
ot 



d(p\) 
dt 



+ V ■ [pVV + p,o,I - BB/po] = , 



(C.l) 



de I BB\ 

— + V- eV + /?,„,V-V-— =0, 

dt \ po 

pV^ P 
— + + 



2 2p„ y - 1 

^ + V-(VB-BV) = 

at 

The density p, velocity V, total energy e and magnetic field B 
are set by the initial conditions as a ID spherically symetric 
SNR blast-wave described by Truelove & McKee (1999). We 
assumed a uniform SNR and added a small contribution of the 
magnetic field. The resulting SNR MHD simulation starts for 
(Ve, = 0) with the parameters 



P = 



3MsNR/poo47rygj^j,r|j^jj 



Vr 



1 

J ^/VsnrT'snr 




,R < VsnrT^snr 

,R > VsNR^SNR 



,R < VsnrT^snr 
,R > VsnrT'snr 



For each run, the physical quantities entering the problem are 
normalised by the known mass ejected Msnr, the age of the 
SNR TsNR, the mechanical energy of the explosion Einj, and the 
velocity of the blast wave Vnr. We set the thermal pressure to a 
small value compared to the kinetic energy of the SNR (typically 
lO"-' times), since its role in the wave propagation is minimal. 
The magnetic field advected along the flow is also believed to be 
very ineffective in the wave propagation but its role in the supra- 
thermal particles transport process is important. The magnetic 
field is thus prescribed with an amplitude similar to its warm in- 
terstellar medium value, e.g.. Be ^ 5jjG. 

To test the ability of our simulation to model the propaga- 
tion of SNR shock, we simulated the long-term evolution of 
a SNR blast wave corresponding to the previous initial set- 
up where we defined the SNR parameter to Msnr = 6Mo, 
TsNR = 200yr, £, 



inj 



10-' 'ergs and Vsnr = 5000km/s. The 
results were found to reproduce the corresponding analytical so- 
lution in lTruelove & McKe3 (Il999h quite accurately. In particu- 
lar, both the free expansion and Sedov self-similar regimes were 
obtained, the transition regime occurring at the expected Sedov 
time for this simulation of Tsedov = 1.1 kyr. 

C.2. Kinetic approach 

The transport of relativistic particles (with velocities much larger 
than the fluid speed) near the shock front is governed by a 
Fokker-Planck equation when these particles resonate with the 
turbulence and enter a diffusion regime. The related kinetic equa- 



tion is 




dF 


d 


dt ~ 


dR 




d 




dp 




52 



dDR 



2Dr 
R 



(C.2) 
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where F - R^p^f is related to the distribution function / in 
terms of the spherical radius R and particle momentum pc - 
ynieC^. The particle spatial diffusion regime is characterised by a 
diffusion coefficient Dr that depends on the turbulence spectrum. 
The factor aioss stands for particle losses. 

For electrons, the losses are produced by synchrotron cooling. 
The cooling timescale tsy„ is 



1 



6nmic 



2^2 



•*syn 



(C.3) 



For protons (or ions), the losses are produced by the gener- 
ation of magnetic fluctuations and are a priori limited to the 
upstream medium (in the downstream flow the particle dis- 
tribution is isot r opic). The cooling timescale is adapted from 
IMarcowith eTall (l2006l) their Eq. 1 3 



'Jturb 



P{P) 



(C.4) 



where P{p) is the rate of energy radiated by a relativistic particle 
Pip) 



1 glog(/(x)) 



(C.5) 



The scattering centre velocity is close to the local Alfven veloc- 
ity, i.e., ^ Vau- 

Stochastic particle acceleration is represented by the energy dif- 
fusion coefficient Dpp - V^p^ /9D^, which is related to spatial 
diffusion (where Va is the local Alfven velocity). 

C.2.1 . Stochastic differential equations 

As shown by iKriills & Achterber3 11994), this Fokker-Planck 
equation is equivalent to a set of two SDEs that can be written as 

dWR 



dR ^ dPR 2Dr 

dt dR ^ R 



dt 



1 dp^Dj, 
p^ dp 



dWp 



dt 



4 



2D, 



pp 



where W, are Wiener processes for which dWi oc -sfdt. Using 
Monte Carlo methods, it is then possible to time-integrate the 
trajectories of a sample of test particles in phase space and to 
reconstruct this distribution function, provided that the number 
of test particles is sufficiently high. 

A shock discontinuity may lead, according to the MHD Rankine- 
Hugoniot conservation laws, to a discontinuous magnetic field at 
the shock front. Depending on the diffusion regime affecting rel- 
ativistic particles, this may lead to discontinuous diffusion co- 
efficients that can be written Dr = Dr c + ADRsign(7? - R^h), 
where the first term is a continuous function. In this case, the 
usual Euler schemes are no longer valid, in contrast to the stud - 
ies of iKriiUs & Achterberg (1994'), 'Casse & Marcowith (2003*), 
Ivan de r Swaluw & Achterberg (2004), and Casse & Marcowith 
(l2005h . As shown bylZhang (2000|), it is possible to overcome 
this problem by employing a skew Brownian motion where an 
asymmetric shock crossing probability is considered. In this 
framework, the spatial stochastic equation becomes 



dR = m{\VR + 



dDv c \ I 



(C.6) 



where R is related to R by 



R = ^(R)R with ^(R) 



R<Rsh 

R = ^sh 
R > Rsh 



and where s is the ratio of diffusion coefficients at the shock 
front, namely 

DA) + Dd(/?sh) ■ 

Equation iC.6\ can be solved using an Euler scheme where the 
stochastic variable Wr is c omputed with Monte Carlo methods. 
In contrast to the study of IZhand (l2000h . reaUstic diffusion co- 
efficients are likely to depend on particle energy. In this case, 
we have to consider the amount of energy Ae gained by parti- 
cles during the shock crossing. The transition probability s is 
then calculated depending on the way that the shock is crossed, 
namely 



^up^down 
^down-^up 



DJRsh, e) 



£)u(^sh,e) + £>d(-Rsh,f + Ae) 
D,(R,h,e + A£) 

Du(/;,h,e + Ae) + Dd(/;sh,f) 



(C.8) 



We note that this skew Brownian motion approach is valid only 
if shock curvature terms are negligible, i.e., 2Dr/R <k IVr + 
dDn/dR]. In the energy stochastic equation, the velocity dis- 
continuit y can be numerically tre ated using an implicit Ricatti 
scheme ( Marcowith & Kir5ll999i) . Basically, once the stochas- 
tic displacement AR is calculated, we can calculate the energy 
gained Ae by a particle originally of energy e - pc during time 
step At following 



e-l- Ae 



exp 



\^ 3A« JR 



R+AR 



V -YdR 



) 



1 -t- e exp 



aiossdR 
(C.9) 

The previous implicit calculation is valid for any diffusion 
regime provided that second-order Fermi acceleration is negli- 
gible. In the opposite case, we then have to step back into an 
explicit sc heme to take i nto account the skew Brownian motion. 
Following IZhangl(l2000h . the energy gained by a particle is 



Ae= ^ 



2D^„dW„ - 



AV 



e{AR-AR/^{R)} + 



dD, 



pp 



- flio 



Af , 



3ADr^' " " \ dp 

(C.IO) 

where AV = yup(-R,sh)- Vdown(-Rsh) andAD^ = DA)-Di{R,^). 
During the time integration of MHD equations, the SNR shock 
front propagates in such a way that its surface increases with 
time. To take into account the increase in the particle flux at the 
shock front, we continuously inject new particles of energy einj, 
so that the number of new particles is A^part(f + Af) - A^part(0 °^ 
R^^(t)ARs\i, where ARsh is the shock front displacement occur- 
ring during Af. 

C.2.2. Kinetic description of MHD sliock waves 

The SDE formalism is useful for modelling the transport of rel- 
ativistic test particles a non-relativistic background fluid since it 
provides both the spatial and energetic distribution of particles. 
Nevertheless one drawback of this method does exist: the shock 
thickness. The SDE algorithm is based on the use of fluid ve- 
locity divergence to mimic particle acceleration. The MHD code 
provides the velocity field at discrete locations on the grid so 
that V ■ V may be obtained by means of linear interpolation. The 
most efficient MHD code cannot reproduce shocks as sharp dis- 
continuities but rather displays velocity and density variations 
over two or three cells. This is very important for kinetic com- 
putations since particles with diffusion coefficients for which the 
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diffusive step is small compared to the MHD shock thickness 
will see the shock as an adiabatic compression, leading to softer 
energy spectrum. 

In previous work (see e.g., iKriills & Achterber^ (Il994l) : 
ICasse & MarcowithI (l2003h ). it was shown that the SDE formal- 
ism was able to describe accurately the transport of particles 
with diffusion coefficients greater than AXshV/2, where AX^^ is 
typically the cell size in the MHD code. This constraint greatly 
reduced the range of applications of this method. To overcome 
this problem, we designed a SDE algorithm in which the V ■ V 
is no longer calculated locally but we instead integrate the term 
V ■ YdR in Eq.( |C.9t , where the velocity is given as Vu or Vd de- 
pending on the shock position. In our new approach, the MHD 
code now provides the shock position and the compression ra- 
tio r so that we deduce the value of the fluid by considering the 
shock as infinitely thin. 
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